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.  Part  i;  inlrocluctibri 

The  possible  scenario  of  clevelopiuent of  chcini<"il  cxplp^ioii  Jnsifle  ciosgd  vessel  resultiuff  in 
a  breakuij  of  the  vessel  Were  described  in  Pn'liinin^y  h’.port  Xo.L  It  Was  shown*  there  that  the 
type  opbrealcup,  the  iiiunlsbi'  and  velocities  of  fragnients  are  detofinined  l»y  tlie  rates-of  energy 
release  mid  wall  pressure  evolutioiiJnsidethe  vessel.  The  re.siilts  of .bfes^iluip  -  mass  and  velocity 
distributions  of  fragnients  -  Wquld  differ  gnxitly  dep«nulihg  (,n4he  cfuhbustion- process  in.sid'e  the 
ve#ei:  deflagfatiqii,  detonatioii  or  deifiagfatioii  to.  diiU>nationuransitiph  i  DDTT 

To  obtain  a  cldsecbforni  splutipn  of  the  problem,  that  ‘’an  lid  .compared  with  the  e.xisting 
expenihent’siweil  try  tP  coiicentrate  oiir  effortsoii  .the  solution  foi  simple  geoiuetrics:  spherical 
^d  cyiihdricai vessels.  The  inithil  coilditipiis  willbe  also  introduced  iii- a  fofnwnaldhg  it  possible 
to  obtain  splutipiis-  ppsassihg  ;.s'piienpahpf  cyhiidncal-s yihhietf y. 

The  results  of  deflagfatipn  to  (letonatiph'tfunsiti()u  investigations  (l.lj  shou'  that  the  process 
is -complicated  and  the  highest  pressures  and  intensities  of  detonation  wave.s  tire  reached  in  the 
trmisitiph  zone.  This  zone  is  characterized  hy  the  existence  of  overdriven*  detonation  waves  that 
slow  clown  gradually  to  self-sustahiing  regime  mul  the  intensity  ()f  the  Waves  decrease  to  that  of 
hbii'iial’  detojiatibh:  ( close  to  Cluipthaii=:T()ug(<t-.fegiiud )* 

hi  case  of  central  ignition  of  spherical  or  ('ylindficul  volume  (‘bnibustiiile.  gas  there  may 
pccur  different  reginies  that  can  produce  different  Wall.-loadiiigs  after  the  refiecti<m.  If  the  DDT 
process  does  hot  take  place  inside  the  vessel  the  pre.ssure  riz('  and  the  wall  loading  is  monotonous 
and  it  is  determined  by  the  rate of  comlmstion:  If  the  DDT  process  takes  place  inside  the  vessel 
rather  far  from  the  wall  or  there  is  a  ilirect  initiation  of  detonation  tiie  loading  of  the  wall  is 
determined  by  a  reflected  normal  detonation  wvve.  Biit  in  case,  the  DDT  ■  |)roces.s  takes  place 
hear  the  wail  the  pressure  rize  (Wall  loadiiig)  is  determined  by  the  reflection  of  an  overdriven 
detonation  wave  that  produces  higher  p(*aks  of  pressure.  It  is  necessary  to  produce  a  closed  form 
solution  of  the  problem  of  wall  lojiclihg,  describing  the  multiplicity  of  scenario  of  DDT-pro(resses 
that  wouldn’t  be  sophisticated  and  could  be'pbtaiiicd  within  rea.souable  proce.ssor  times.  To 
solve  the  problem  it  is  necessary  in  addition  to  the  DDT  investigations  (1.1)  that  match  the 
fe^tUs  of  the  experiments  in  tubes  (1.2:  l,3j  to  examine  some,  peculiarities  of  spherical  and 
cylindrical  deflagration  and  aeloiu.tiou  waws  propngntion.  Part  2  c'ontains  the  re.sults  of  these 
ihvesligatiohsi 

The  clo.sed-form  inocWl  of  wall  loading  basc^d  on  the*  ic*.sult.s  of  investigations  conducted  taking 
into  acciount  the  multiplicity  of  scenafio  of  DDT-procc'sses  is  described  iii  Part  3. 

Part  4  describes  the  results  of  Ijieakup  modelling  tmder  the  influence  of  internal  loading  for 
different  initial  conditions  and  different  sceusirio  of  comimstion  process  inside  the  vc'ssel.  It  is 
shown  that  the  criterion  of  breakup  ctin  be  determined  from  an  iiulependent  experimental  data. 

The  mathematical  mcjdel  of  breakup  takes  into  ac-count  the  influencc.'i  of  temperature,  on  the 
behaviour  of  the  main  i)araihe.tefs.  It  is  sliown.  that  tlu*  critical  pressure  of  breakup  changes 
depending  on  the  type  of  combustion:  the  higher  is  the  rate'  of  loading  -  the  higher  turns  to  he 
the'critical  breakup  pressure.  The  number  of  frtigmcmts  is  higher  for  higher  pressures. 

Part  5  presents  the  results  of  numerical  modelling  of  fragments'  acceleration  in  the  expanding 
cloud  of  gaseous  reaction  prodnerts  aften-  the'1)reakup.  The  immericcd  c'xperiment  conducted 
made  it  possible  to  Work  out  tipproxinuitiou  formulae  di'termining  flnal  velocity  of  fragment 
depending  on  its  .shape,  orientation,  initial  conditioirs  and  thcninodynamic  properties  of  the 
expanding  gas. 
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wayes,  with  spherical  syhl- 


^  At.deflagratipivqi:  detqnatiou;  imtiatiqn.iu,  preiuixc'cl  hoinogeiieous-  (■oiubustible  liiixtiues,  for 
.  exaiyiie  by  .meaiis  ()f  a  siKirkpr  a  detonator.  a4:eiativ4y  .siiiaU  zone  of  raised  values.of  parameters 
. appears*. t^sizpne  caii  . be, 5ionsi(iei^^^  a  point  s(mrco.at  a  .siifficient  distance  front  the  hiitiation 
.  location:,,  .Expan^^^^  proclijcts,  ,of  detonator  e-xplosion  In  .sonie  (’.a-se-s  cap  play  the  nile  of  a 
r- splieric^:  ^  •  \ 

Self^sinijiaiity  of  the  .pipblein  oji  sphericjii  dchtnation  wqs  estabiished  and  the- ordinary  dif- 
ferentiai;pquatipnbv.f^^^^^  describing,  the  How  behincUa  wave  front  were  tterived  in.  j'il).  The  case 
,qf  tlie„  ^haptUanr  Jouget  sphericai.  detonatioji  \\'a.s  (tonsi(l«r(’d  hir  exainple  in  (2.2.  2.3).  The 
fro§t,cpmp|$.te  splutip,nrpfr|he  .self-simiUw  proliiei^  case.s  of  .wetik  and  siro  ’g  fletonation  and 
.Ghapmtui-lpuget.  tyayes  was  oi)taiiied  an  {2;4j.  In  the  present  part,  j)roHie.s  of  How  parameters 
m  propagation  of  .sphericai  w<aves  <>f  strong  aAiuiwefdv  dettmatioh  and,  deflagration  for  cases  of 
different:  yeio,bitie^  of  .iHsttm  nujtio^  from  the  initiation  zone  are  inv!e.stigate.d  on  the  base  of 
imethodsAdlyelpped’^.m^^ 

\Ye  uge  . the  model  where  a  detonathni  ancl.dbflagfaticm  >vave.s  zone  i.s  asstimed  to  be  a  discon¬ 
tinuity  surftice,  and  suppose.aiiat  detonation  arises  at  tihie  t  =  ddii,  the  symnietry  center  and 
propagates  with  a  coirstaht  velocity:  the  initial  ani.xtiire  and  reaction  products  are  polytrbpic 
ga^s;.^o.#cohst  i.s  the  specific  heat  of  reaction;  C'ohil)nstil)l(*  ini.xttire  is  iionniobile  ahead  of 
the  waye,  ;and'entfppy  of  ali  .particles  has  a  constant  valiu'.  and  entropy  incnnuses  for  all  parti- 
byfhe  same  yahte.  in  the  detonation  wave- moving  with  constant  velocity  U  i.e.  a  How  of 
detonation  products  within  the  continuity  domain  i.s  lipnioentropic,  Etiuaiions  of  nonstatiohary 
on^diniensibhal  motion  ha.s  the.  f  dlowiijg  Hnan  under  the.se  conditions: 


dit  Ou  ^  ..i  Qli, 
Ot  *[dv  i)0v' 


dp  Op  Ou  'Ipu 


=  0: 


(2.1) 


where  a,  p,  p-  ^iis  velocity,  pressure  and  (h'usity  re.spectively.  Under  these  a.ssumptions  eciuations 
and  boundary  conditions  allow  an  existence  of  .self-similar  .solutions.  i;e.  solutions  being  functions 
only  of  one.  parameter  (  =  r/f.  P.ussing  to  self-similar  varialde  f  by  the  formulae: 


Ot  t  (/f  ’  Ov  t  ’ 
we  reduce  equations  (1)^  (2)  to  the.  form: 


(III  \dp 

—(a  -0  =  ---77' 
f/f  p 

(2.4) 

{ti-0<h>  .  .  2(f  „ 

(2.5) 

Since  the  flow  is  homoentropic  then  dp/d^,  —  c-dpjdi,  where  c  i 
equation  (5)  to  the  form: 

s  the  sound  velocity.  Transfonii 

\dp  idu  ^  '2a  ^  (•' 

(2.6) 
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should  i^e  nptifeetl^  a  coustaut  flow  Ji,(/.=c()ust|  isirt  a  solutiou-of  equatipu  (2.X)  for  the 
Ijeneral.c^e,  Tri^flal  scflution  of  ii  (:p.its.u  i.e,  u  =/(],- (.a  quicsf'ent.  state)  is  ouly  i>ossil)le; 

Let  us-seek  a  iio,n4riyial;seif-si,niii;ff  spiulitni.  linrofluce  iiew  (h'peiqlent  yaiiaiiles 


:andtihcleii5e;ucleiit  viiiiubW 


'/ ;='lh,(>^;). 


gping.'.tp' ihe  new  vaiiai)le.s  the  system  l’i.G);  (2.7)  rakes  the  foi'iii: 


where  F  =  For  poly^ojtie  gtis  F  =  (-  -  1);  where  *)  =  the  polytropic 

exphiieiit.  Itfisplving  (S),  (9)  for  the4erivative.s  we  obtiiin; 


(hr 


(2.10) 


(2.11) 


With  the  lielp  of  the.  .system  (2J()).  (2.11)  we  derermiiK'  ^  and  then  write  the 

expression  derived  in  place  of  etpiaticm  (2.10).  Thtts.  the  following  system  is  derived: 


ch  2::c-(r-.l)(yrtA) 

(IV  ^  V  3c-  (l-^F)’i  ■ 


(2.12) 


fh  1.  ~-(l.^V)'  /.)13) 

(H'  F3;-(l-r)-^' 

Sy.stem  (2.12),  (2.13)  can  be.  solved  setmentially:  tu  first  equation  (12)  gives  the  solution 
z  =  s(y)  and  then  the  solution  obtained  is  substitutetl  into  (13).  from  where  r;  =  (jiV)  can  be 
found. 

Define  boundary  condititms  for  .sulyurg  the  etpurtion  ( 2.12)  in  the  phine  ( V,:).  For  fixed  t  the 
condition  ^  =  ^(r)  is  satisfied.  Coasider  the.  zone  ahead  of  the  detonation  wave.  Bo  <  r  <  .I 
(Figure  2.1).  In  this  zone  V  =  ii®/^  =  0  (a  quie.seent  zone).  In  this  (;a.se  ^  changes  from 
=  oo  for  r  -*  oo  to  —  R/t  =  U  for  »■  =  /?  =  Ut.  Value  of  z  within  this  zone 
changes  from  zero  (c,.i  =  0)  for  ^  »  oo  to  h,  =  cijU’  for  E,Bo  =  On  the  plane  (V,.:)  it 

corresponds  to  a  trivial  solution*  when  a  state  changes  from  iroint  -4  (F  =  0,  =  0)  to  point 

Bo  {V  =  0,  c„  =  4/Lr^  <  1 )  (Figuir  2.2). 
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.  ipsttmnetersJjehiuU  thtwie  xyiiyp  tiphii  B \  ) » an  be  (Ictfi-miiKMl  lioiu  the  (equations  at 

fch(ii^iscQntmui^^  .  ■  ',  '  ,s? 


,  •  ■'? 


•?> 


. .  r(Ui  ^  i?0' T /'()  f' n  • 

.  ..r.'  ,,  --ei.  ,  ■r',  .  ,  / 

—'T.  i==;0;i': 


tfi^e  cduciUkhi(>4pv  t.lie.ia'ojjiniu  .nnisi  at  ll<  =  0  tak(‘,s  t.h(?  ,fonu  [2.4]: 

.  ■  Pi 

-I .+  ^  +iiA); 

;u 


,  *) 
Six . 


where 


*■  /I  T  1  ;«  J 

v-i  - 1  _ 

'  7ra  +  ^) 


(2.16) 


(2.1S) 


It  fdliow.s  from  (2;i6,)'2mtl  (2.17)  that  for  fixod  Aft')  tlu*  point  D\  in  the  phuie  (F, r)  is  oii  the 
ptuabola 

1  _L  -»,i\ 

(2.19) 


4,  = 


l^A 

At  the  Chapman- Jonf^et  detouatiou  it  rpiaiiye  ga,s  t-oloeity  Inihiud  the  wave  is  equal  to  the 
local  sound  velocity: 


Thus,  tlxe  Chapman-Jouget  detonation  eoiTe.spond.s  to  the  case  =  0  wlu'u  the  parabola  (2.19) 
takes  the  form 

c,  =(l-r,)-  (2.19(0 

and  goes  through  the  points  (0,1)  and  (1,  0)  (Figure  2.2). 

For  >  0,  U  >  Up  tstpuition  (2. IS)  has  two  roots:  .V^  >  0.  A"  <  0;  |A“|  =  A"*".  It  follow.s 
from  equations  (2.17)  and  (2.15)  that  >  nj',  p'l  >  />]*  i.i*.  the  root  >  0  coiTe.sponds  to 
the  case  of  strong  detonation,  and  the  root  A"  <  0  corresponds  to  the  case  of  weak  detonation. 
Solution  on  the  plane  (l<^  z)  can  be  prolonged  (dther  to  the  z  axis  (i.e.  V  =  0)  that  corresponds 
to  the  rest  in  the  symmetry  center  dr.  at  pivseiice  of  a  sphtnical  piston,  expanding  with  rate  t(,„ 
up  to  the  line  F  =  1  as 


»•;.  =  It  J 


t 


f 

I 


In  the  zone  between  the  piston  and  the  (U^*tonation  wtive  ti  shock  wave  can  appear  for  which 
the  same  relations  at  a  <liscontinuity  stnface  (2.14)  are  true  tus  for  ti  detonation  wave  but  with 
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cphdilidii's  that  there  is  no  eu^rsy  release  ilny  r(),a  cheiuicai,  leactiun.  i.Qo  =t  *))  •‘nh  ‘V  luixnire 
hpiiiiypsitiou/thies^  irhan»e  (  e:^  •  l(T:)  ;=  t  T j); 


ihiM  /'iC!  j  i;)? 


,  Til ., 


hi  —  i  +  'T  1  V  —  ^*1  ~  ^  -  “(I"-  -  iv' 

mV-y^l)  7i(hi'-l) 


;V^.==  1'^ 


:2V,, 


=  m  -  i)" 


..:  .  ,1.-1-  1,  =  i 

Gbnsider  ihe  case  (h  ,Chapmau-.|(mg(H  spherhral  (l<\t<)iiati()ir.\vav('.  pnipasatioirlX^  =  Uq,  -^  = 
-0.).  Sinje  ^v,.>  po,  it  foilows^rQui  wiuatuni  (,iip>  tliir.t  U  <  Vy  <  1  -  jH)Ipi.  <  1  aiul  that  the 
p.pint  B,i;  (y  =  Ki.,  -  =  (('.\  /Up)'  <  1),  is,,()u,the  left  l)raueh  of  parabola; (  2;i9,a).,  ,Pa.ssage  fnyu 
point  J?o  to, pbhit  B\  in  the  plane  {V,  c) ocaan-s  by  a  jump  (Figure  2.2);  In  the  ca.se  wlieii  gas 
in.-the.symmetiy  center  is  at.,re.st,,(M,  =  5;  r  ='0)  point  M  (I^igure  2.1)  in.the  plane  (y  i)  is 
on  the  twis.y  =  0:  M(^Vsp=  d,.2,\/  ?=  co);  Point  jv{VV  =  O.c/^-  =  1)  is  a  .singular  point  of  the 
knotstypeifor  equation  (2.12).  Thus;  a  .solution-cOnneetiug  point  B\  '.ith  pohit  M  goes  through 
.point  . 'A';  Th’isjpoint  fiKin-niaterial)  UMives  on. the  iilane  (;•./)  wuli  the  sound  velocity: 


“t-  =  ^  =  1 


's- 


n- 


r  =  ryf: 


n  = 


iii 


—  0  (ti  =  0. 


i.e.  there  is  a  characteri.stic  on  the  plane  (i\t)  which  boinuls  the.  re.st  domain  (at  this  charac¬ 
teristic  derivative.s  have  discontinuity).  Solution  Irom  point  A'  to  iioint  .V/  is  re.st:  as  in  point 
A'  the  condition:  Uk  =  0  is  satisfiwl  and  in  .point  ;)/  i/;i/  =  0  dpe  to  symmetry,  then  the  rest 
.state  satisfies  the  equatiims  and  the  boundary  couditibns  and  due  to  uniqueness  it  is  a  solu¬ 
tion  between  the  points  K  and  i\A  On  the  jilanc  (V'c)  states  from  K{Vi^-  =  O.j/y-  =  1)  to 
M{Vm  =  0,  =  oo)  W’**'  the  c  axis  from  1  iip‘t<>  .xj  (Figure  2.2),  Thus,  a  de.sirecl  solution 

behind  the  detonation  wave  consist.s  of  two  Itranclies:  Inanch  KM  corresponds  to  the  trivial 
solution  y  =  0,  and  branch  A'G’oiJi  to  the  .solution  of  the  general  type  for  .sy.stem  (2.12),  (2.13). 

Oh  the  plane  (r.t)  a  .solution  bdiind  the  d<‘ionatiou  wave  (/■  =  L'ot)  consists  of  the  expanding 
domain  of  gas  at  re.st  (a  =  O.p  =  con.st  =,/)^j  and  the  rarefaction  wave  in  the  domain  b(*tween 
the  characteristic  C'|':  v  ^  and  the  detonation  wav<*  r  =  Upt-  This  rarefaction  wave  isn't 
a  solution  of  a  simple  wave.  .‘Moug  each  axis  ^  =  covst  values  of  How  jjarameters  are  constant. 

Consider  the  ca,se  of  strong  detonation  (.\  =  A+  >  0,  U  =  U\  >  Uph  "'hen  gas  goes  from 
the  state  characterized  on  the  (V,z)  jilan  ■  by  point  i?oi  to  the*  state  characterized  by  point  Bn 
(Figure  2.2). 

In  this  case  point  Bn  characterizing  pjirameters  immedi.-itely  behiml  the  detonation  wave  is 
situated  oh  parabola  (2.19)  which  is  higher  then  parabola.  (2.19a).  Coordinates  of  point  Bn: 
V  =  Vii  and  z  =  can  be  determined  froin  eipiations  (2.1G),  (2.17)  (Figure  2.2).  In  this  case 
we  can  construct  the  solution  where  a  comitrrssion  wave  propagates  behind  ti  detonation  wave 
(it  is  an  analog  of  a  constant  How  in  a  phiue  ca.ser(2.5l. 

Continuous  solution  B\\Pi  (Figure  2.2)  n'aches  the  line  V’  =  1.  points  of  which  can  lie 


considered  as  points  corn' .ponding  to  spherical  pistons  expanding  with  constant  rate  Point 


Pi,  wherein  the  integral  curve  comes,  defines  the  only  \'elocity  of  the  piston  which  corresponds  to 
the  given  velocity  of  .strong  detonation  wave  B,.  With  decreasing  Bi  parameter  A"*"  diminishes, 
and  A"^  =  0  at  Bi  =  Bp,  parabola  (2.19)  transfers  into  parabola  (2.19a)  and  point  Bu  becomes 
coincident  with  point  B\.  Continuous  compre.ssion  wave  is  described  by  solution  B\Pd,  the 
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4MStph  yeldcity ^decreases)  Ip  value  i.e.  \ye  obtain  a  solution  with  the  C'hapiuaii-JousLet 
detonation  wave  aiid; the  oontiiuipuscompressionr  wave  liehind  the  detonatioii  wiive  front  to  the 
iiphmctd  pistph  nioyiii^  udtlV  velbiuiy'M^^^^  , 


On  reducing  the  pistoii,  velocity' so  thu’t  a  state  ahead  of  the  piston  wciuld  <’orrespond  to  point 
1P2;  in  Figure  2.2,  we  get  the  solutittn  containing  the  Clitipinan-Jouget  detoiiatioii  wave,  the 
rtnefactiph  wave  traitsfPrniing  ga.s  froiii  state.  By  into  certain  state  C\,  on,  curve  B.\  C()K,  the 
.sliodi  wave  tran.sforining  gas  front  state  Co  into  state  C'l  oil  curve  B\C\K.  and  the  continuous 
cpinpressioh  wave  C\B>.  Curve  BiC\K  presents  a  locus  of  states  on  the  (!'  plane  which  caii 
be  obtaiiidi  as  a  result  of  a  shock  fiausitipii  front  the  states  lying  oil  iiitegrai  curve  B\C\jK. 

With,  futiier  tlecreasing  th<*  pist<»ii  v<‘locity  u,,  coprdiiiate  r,,  of  poiiit  P  grows  up.  In  this 
ca^  there  is  a. .solution  wherein  points  C\  aiid  C'o  iiiove  along  curves  iJiG'ij A'  and  B\G\K  aiid 
approach  the  point  K.  Fpf  u,,  =  ()  there  is  a  .solution  containing  CliaphiaiirJouget  detonation 
wave,,rMefactipii  \vaye  aiid  rest  Zf)ne.  Under  the  giyen  houndary  (’oiiditionsf  at  the  pi.ston  u  =  n,, 
and  at  infinity)  the  constructed  solution  hi  the  cltiss  of  strong  detoiiatioii  aiid  Cliapinan-Jouget 
detphatiph .  ( A  >  iO ) dsn ihique; 

Figure  2.'3  shows  computed  fe.suit.s  for  the  case  of  spherical detonation  iii  mixture  Qo  =#  3  •10" 
.I/kg;  7o  =  U4;  71  =  1.2  in  the  form  of  *de.jk*iid(nic(*  of  pressure  /i/po  on  noh-dimensioiiai 
cooidinate  r/R.  Curve  BtyBfCoKM  cofre.s-iionds  to  the  solution  I'Oiitaiiiiug  the  Chapniau- 
.iouget  detonation  wave  B(,Bi  (Un  =  1030  ni/s),  the  rarefaction  wave  and;  tile  rest  zone  at 
itj,  =  0.  Curve  Bt)13\C\tC\P,  eorfesponds  to  tin*  .solution  eoutainiug  the  Chapiiiaii-.Iouget  wave, 
the  rarefaction  wave,  shock  wave  CuC\  and  compre.‘<siou  wave  C\P>  at  the  expanding  spherical 
piston  velocity  u,,  =  1200  iii/s;  Curve  B(tB\CliG\P',  is  the  pressure  profile  at  decreasing  the 
piston  velocity  itj,  down  to  1000  in/a.  It  is  .seen  from  the  figvire  that  futher  decreasing  the  piston 
velocity  will  lead  to  reducing  an  intensity  of  .secondtiry  wave  C'uC'i-  Curve  Bi)B\Pi)  corresponds 
to  the  .solution  containing  the  Cliqmiaii-Jouget  detonation  wave  and  compres.sion  wave  B\Po 
at  piston  velocity  w,,  ~  1320  m/s. 

Let  us  go  to  consideration  of  weak  detonation  ( A”  <  0).  In  tliis  ca.se  tlie  prolilem  solution  isn’t 
unique.  An  unique  solution  can  be  olitauied.in  the  weak  detonation  class  wlien  a  value  of  one 
parameter  on  tlie  detonation  front  is  given  in  addition  to  tlie  boundary  conditions  (2.4,  2.5),  for 
example  this  pai’ilmeter  can  lie  detonation  wave  velocity  t’j  which  define.s  uniquely  A7  <.  0  and 
point  Bu  on  the  {V,z)  plane.  Solution  contaiiis  the  rarefm'tion  wave  By>Co'>K  (Figure  2.2)  aial 
the  expanding  rest  'zone  liehind  the  weak  di'tonatiou  at  11,,  =  0.  Pressure. profile  for  the  solution 
obtained  at  17)  =  1750  m/s  is  shown  in  Fig.2.3  by  curve  BuB\  j\M.  Solution  at  presence  of 
a  moving  piston  («,,  =  const  >  0)  contains  the  rarifaction  wave  behind  the  detonation  front 
from  5i2  to  certain  point  Cu*  lying  on  integral  curve  By>G{yiI\.  the  shock  wave  transforming 
the  medium  from  state.  Coa  into  state  Cy>  on  curve  B\\CyiK  ami  the  continuous  coiiipressiou 
wave  to  state  P  {V  =  \  )  ahead  of  the  piston  (Figure  2.2).  Curve  B\ \  Cy>K  in  Figure.  2  is  a  locus 
of  .states  which  can  bvs  reached  by  a  shock  transition  from  slates  on  curve  BnCo^A'. 

Figure  2.4  shows  pressure  profiles  for  ih(‘  case  of  spherical  weak  detonatioii  wave  propagation 
with  velocity  U->  =  1750  m/s.  Curve  B»\By>CMiKM  corresponds  to  the  ca.se  xip  —  0  when 
there  exists  a  rest  zone  in  the  symmetry  center.  Solutions  with  an  expanding  spherical  piston 
are  represented  by  curves  B{}\By>Cw>Cy>P  {up  =  1300  m/s)  and  BinBooGyytC'^oP'  («;>  -  1000 
m/s).  As  seen  from  tliis  figure,  witli  increasing  the  piston  velocity  tt,,  an  intensity  of  secondary 
shock  wave  C02C12  grows. 

Solution  of  the  indicated  type  for  the  ca  .e  <  0  exists  not  for  all  values  of  the  boundary 
conditions  but  only  till  the  piston  velocity  i/j,  <  where  u,,\(U>)  is  the  piston  velocity 

satisfying  the  soliition  with  a  strong  detonation  wtivi'  moving  witli  velocity  U\  =  Uo.  Curve 
BoiBnBiiPi  in  Figure  2.4  correspomls  to  a  pre.ssure  profile  for  such  limiting  case  (  For  tlie  case 
considered  =  «;.i(£72)  =1500  m/s).  .\t  this  piston  velocity  a  secondary  shock  wave  C'o2Ci> 


qyertc^es>A*#rong,  (ler<)iicvti^^^^  uwi;gqs  witli  that  anti  tbiins  a  strong  detonation 

waye,Bor5i!i;.nhqying  with . tin?' sanie  yebcity  f'i.  =  1100  in/s.  For  gieater  velocity of  the  piston 
Up.  >  ,Ujp|C/r2)^a  stVkitio  T).e  fpinnl  only  hi/the  class  of  stipng  detcin^tion,  For  tliis  case  a 
;iiecess|iy  to  giye^an  adcjitiqnai  paraihe.ti'i-  at  a  de.ton.ition  iroii't  is  eiiniiiicitecl  (2.5). 

Joilttioxis  of  atiaiqgic5ii-,prtd/k  detonatiph  \v:iv(?.s  motion  with  the  additional  effect  of  a 

,pistqh*haVe  beeiv-olttainedrftn'  the  jdinie  piiertlim^^  case,  for  i?xaiupltj  iii  f2.5j. 

Gdnsider  the^pnihieni  tin  .spherictil  Hami*  front  propagiifion  from  a  point  source  of  initiatioii 
."for  the ‘case  ydieii  deflagratmn  (hxesiFl  in.'conu;  .dtuonati.om  DeHagratipn,  front  is  assumed  to 
iief.a  discpiitinlutj'  siuface;  aiul  t  {c.s,()lutioh  is  being  soiiglit  within  the  scope  of  the  siiine 
^ais^iuptipus  as  in  the.  previous  problem.  Coirsider  tlie  cttst'- when.. a  \veak  deflagration  frcnit 
^iropagates  with  a^given  constiuit  rate  tfl'  inunud  buniing  ii’  throngh-.au  initial  quiescent  nnxtiife 
AvitLc6nstaht;pai'aineters.  Nprinal  burning  rate  tij' is  siil)s()ni<;  withaespect  to  gas  ahead  of  the 
frpnti  therefore  a  scdutiqn  \yili  contain  a  .spheiicai  shock  wave  moving,  ahead  of  a  flame  front. 
To  finti  a  solution  in  the  disturbed  zoiie  of  nou-comliusted  gas  hctween  a  shock  wave  aiid  a 

deflagratipniront  tile  system-(.2.i2h  (2.431  caii  be  used  tyith.cfnresppm.iing  inniudafy conditions 
onlhe  shoch  waye-  and  onrthe  deflagruticm  fnjiit- which :hav(*  the  followiirg  form  !n  variables  ( K  ^): 

liiOipthe  shock  wave  nioving,thf<mgh.(|uiesceut,nie(litnu  t  \i(  =  0).fr(>m:(2.15)-(2.i7)  we  have: 


=Ti=lq.)(lp 


•r,). 


{2:21} 


' '  I  .  . 

ffotir-where  it  is  .seen. that  slnn^k  wave.s  tninsforni-ilie  tixis  i’  =  0  into  points  of  paraiiola  (2.21): 

2)  oii  the  flanife  front,  behind  which, a  state  reaclies  the  one-  at  rest,  parameters  ahead  of  the 
front  are  connecteddxy  the  relation: 


(2.22) 


(  l.f"  ^  7  “d'V)'- 

70  70-t-'^  •’ 

iiC.  going  onto  the  axis  V  =  0  throtigh  the  flame  front  is  possilde  only  from  points  of  curve 
(2.22)i  in  deriving  the  equation  (2.22)  we  dse  etmditions  (2.10)*  (2.17)  at  the  discontinuity 
.suffjice  and  the  relatioii  connecting  a  normal  burning  rare  11’  and  flame  velocity  C^/;' 
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Parameter  70  corre.sponds  t<>  fh<?  initial  mixture.  71  corresponds  to  the  comlmstiou  product.s. 

The  zone  ahead  of  the  .sho('k  wave  images  into  tlu*  s(«gmeni  AH(,  on  the  (K.  c)  plane,  where 
coordinates  of  point  Hu  corresponding  to  the  .state  :ihead  (>f  the  shock  wave  front  are:  V'o  =  0: 
-0  =  Cq/C/'^  <  1.  The  state  Ixehiud  the  shock  wave  eorre-sixaifls  to  the  point  H]  lying  on, parabola 
s  =  zi{V)  (2.21)  (Figure  2.5).  Going  from  state  H\  into  statt'  :V»  ahead  of  the  flame  front  lying 
on  curve  z  =  (2.22)  is  realized  within  a  continuous  compression  wave  (Figure  2.5)i  and 

traasitioii  into  state  Ni  behind  the  flame  fumt  (Kj  =  0.  ;t  =  ryu}  >  1)  occurs  by  a  jump. 
The  rest  zone  behind  the  flame  front  image.s  into  ray  .Vf.l/.  For  this  case  !i  solution  «.s  unique 
at  weak  deflagration  velocity  11'  given. 

Notice  that  points  of  curve  (2.22),  for  which  coordinati'  ;  $  <  1  at  point  A’),  cannot  corre.spond 
to  the  leading  edge  of  the  flame  front,  as  that  leads  to  s\q)er.sonic  flaim*  front  velocity  with  re.spect 
to  gas  behind  the  front  (strong  fleflagraiionj.  For  this  case  tla*  solution  can  Ije  .  .mstrncied  that 
contains  the  Ghapman-.Iouget  deflagration  wav('  and  tlie  landaciion  wave  lollowing  the  first  one. 
State  beliind  the  deflagration  wave  should  be  situated  on  j)aral)ola  r  =  (I  -  1’)',  and  a  state 
ahead  of  the  wave  is  on  the  curve: 
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presentjing,  a  locus  ()f  states -iyqm,  \viucli  states  ,ou  iKiniljola  i2.ipa)  can  Ijp  iKached oil  goina 
-through  tlie  dame  The  st  put  ion  ohtuii&^cl  contains  the  slioch  transition  frqiii  the  hiitihl 

^ppmt  pn;..a?cis.  V  ,  =  hitc)  . point  H\\.  on  cum?  , then  the.  transition  lyithin  a  cpntiuupus 

Cdmpfessipn,  vwavf  t()  .state.  on  curve  (2.23,),  the  jump  iii  tile  deflagratibh  lyave  to  stSte 
F  on  ptir:abpiaid‘2.lPaJ  fiu^K^^  the.  coiitmuous  transitipu  within  the  rarefacthni  wave  to  pciiiit  A 
(Eigiirp  2,5),.  This, picture  supptises  tv  presence  oha  zone  oTquiesceiit  gas  in  the,  .syininetry  Renter. 
Depending,Pn;aipc.atipn,  of  point  F  pn;ptirabpia  (2.iQa)i  ftirther  mp.tittn  is  possilile  either  aldng 
inte^al^ciiiwe  F&(iA'  tp..p()int  AApr  t<)singi\lai%ppiut-.D  (;  =  V  =  or  to  singular 

point  -For  the  sejunid;  ca.s;e  die  rest  zone  d(><^sn^t  tpjnl;  tiinl  thc  inoiiqn  .ctnitihues  to  the 
synuuetry.ceilter  and  .for  the  last  case  enqjtiness  fqrins  netir-the  cfiiiter.  It  is  nece.s.sary  td  iioyice 
that\there -exists,. a  s()luti(m  fdr  .alhpoints oniparaljqla  (2.21.)  higher  thaii,  point  ‘I  of  pardholtn' 
(2i21).^and;.('2,23),  intersection;  (lilvparaiiqhi  (2,19a)  point  «I»i corresponids  to  this  point  ‘5  and 
the^exmuple  (.ainsicWed  is  Tiigher  than  , ixmit-Fi  ,  that  says  there  are  solutions  vvitli  a  rest 
zdne'fdrAti\is‘Cpe*atvah'senced)f-a;pi^  ,  ,  ,  . 

Sqhitidn  for.  weak  deflagration  (ir  =  3()  in/s,  U  =  500  in/s)  on  the  (i.  E)  plane  is  shown 
in  Eigure,' 2,6.  Segnieu  from  0  to  if/,,  ccariisponds  tt)  a  (juiescent  hiitial  iximimsible  nuxture 
zone;  .Jump  to  point  AT  occurs,  throiigh  the  sluick  wave,  hitegral  ciirve.  ATiV-j  <.’oiTe.sponds  t()  a 
continuous  conipressioii  \vaye  ahead  of  the  flanie  front,  .luinp  of  parivnieters  .y<;V;i  qcctirs  through 
the  deflagration  front.  .Aaest  zone  be.giii.s  (lirectly  hvAind  the  (h'Hagration  wave.  Piessiire  profile 
cdrfespdhding'to  the  scilntioii- obtained  is  .shdwn  liy  eurv(',  ifoAq  A^A';i  in  Figure  2.1. 

Solutidn  with  the  GhapmaurJouget .deflagration  wave  (IF  =  5S0  m/s,  U  =  S60  m/s)  on  the 
(zyV)  piaiie  iii  Figure  2.6  corresponds  to  segment  Oi/ti  (a  ([uieseent  initial  niixtui’e),  jump  H{)L 
tliroughdiie  shock  wave,  c.ompressiou  wivw?  £.D,  jump  thrdugh  deflagration  wave  DF^  continuous 
rarefactidii  wave  FK  and  the  zone  of  quiesceiit  gas.  CorriJspondiug  pressure  profile  is  .shown  in 
Figure  2;7  by  curve.  HoLDFK. 

At  presence  ()f  an  expanding  spliericiil  pistouTor  the  Chapinan-.i(mget deflagration  there  are 
pos.sible  sohitions  with  a  compression  wave  iusteafl  of  a  rarefaction  wave  behind  the  defiagratidh 
Adnt  (integral  curve  FP  ,iu  Figure  2.6)  (ir  with  a  rarefaction  wave,  a  shiick  wave  and  cdmpressifni 
wave  by  analogy  with  the  i.’ousi(lere<l  fiqiy  pictun^s  behind  the  (.’'.uipmai^Jouget  detonation 
front;  Fdfiexample  Figure.  2.7  shows  a,pre.ssure  profih*  for  piston  vt.'ldcity  =  5S0  m/s  (curve 
H<,LDFWM. 

Thus  possible  flow  pictures  at  jiropagatiug  deflagration  iiud  detonation  waves  with  .spherical 
.syihinetfy  have  been  investigated. 


.  ,Pat4^*  of of  sphei’icai.or  cylindrical  vessel 

7pfediffef  of  comfeiistipn  pf passes;  inside  .t  he  vessel 

Qne  of  the  luaiu  paniiueters  of  the  4jfi()rc.ss  ,i.>-ihe  of  ••oinbiistible  mixture  in.'-ide  the 

vess4-  £et.  Its  assiinie  that  the  vess‘el,  routajus  i.'i;:iuixt}U(;  of  i(;a»'.*at.s ;au<l.inob.al)ly  iiiert  compo¬ 
nents.  r^‘  T?mass  concenlrnlioirof  reagents  ( /  =  1  -oxidizer.  /  =  f  -  fiad.  /  =  3  -  iinnt  components. 
/!;=  4  .-d'oactipn4),iodUc.ts);  //,  -  nudaf 'inas.si\s  ()f:tht^  coinpommts. 

The  hfuttorhjactitin  ^can*  l).eidij.scril)ed' l»y  t  h»‘-foiinuJa: 


where.:!/-,  u"  r  suhchioinetiic  coeffiehuits  lief(n%'-m>d.Jifter  tlu*  re/iction.respectiyely; 

If .diemicalenergy  of  e«)thijustioii  p(*f  ma.ss  unit  of.fiiel  is  known  (let  it  he  Ai?)  the  specific 
energy  feiea.se  <5  (eiieiiiicai  energy  per  hiiiss  iluit  of  mixftin''  rail  he  (leteriihn'efl.  ius: follows: 

^  ifr, 

\vhere  <5  =  =  stoidhtjihetfic  ratio-. 

Let  the  total  mass  of  ga.se.s  ..sidt*  a  spherical  ve.sst‘l  he  M.  For  the  cast'  of  cylindrical  ve.ssel 
let  M  bo  the  hiassof  gas  :per  unit  length  of  ti  cylimlef.  The  maximal  etpiililnjum  pressure  p,„ 
in  the  vessel  after  the  comlnistionof  mixnire  can  he  determined  hy  ii  formula  derived  from  the 
energy  conservation  law: 

lii"  =  (3.3) 

Po  ^0  /hi(  '  C.,i  C,,,T1( 

where  T,n  -  avera.gft  temperature  of  mixture nfter.comluistiom  .c^,.  c„  -  .sitedfic  heat  of  mixttire 
before  aird  after  the  renction:  /toi  /On  •  molar  ma.sse.s.of  mixture  before  and  after  the  reaction: 

=  0.4> 

sT''' 

To  obtain  the  clo.se.d  form  motlei  of  walls’  loading  it  is  m-cessary  to  (examine  different  regimes 
of  combustion  proce.ss: 

1.  Deflagration  wave  propagate.s  from  the  centre  of  .symmetry  to  the  walls  (Re.gime  Xo.  1). 

2.  Normal  detonation  wave  reflects  from  the  wall  and  i)rodnce.s  impulse  lotiding  (Regime 


3.  DDT  proce.ss  takes  place  near  the  wall  that  causes  the  reflection  of  unsteady  wave  from 
the  wall.  This  wave  of  unsteady  strueture  can  be  w<‘ak  and  it  can  Ijp  very  stron.g  in  case  of 
overdriven  detonation  (Re.giine  No.3). 

One  of  the  main  characteristics  of  tlie  DDT  provess  is  thi'  pnsletonation  length 

L  =  L{pQ.  T»,  Co.  yij  .  y'j' ....) 

-  the  distance  from  tlie  ignition  point  to  the  phice  of  the  onset  of  detonation.  Experimental  and 
theoretical  investigations  of  DDT  processes  (3.1  -  3.14)  in  closerl  V(\ssel.s  and  unbounded  volumes 
showed  that  the  predetonation  length  depemls  on  many  partimeters:  pre.ssure  (po),  temperature 

I  1  Typeset  by  .4,vi5-'r^X 
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(To’)f'chemicai;'pi()i)eru^^^^^  mixture  kyi'j!'.  j  =  1.2; ..w-  •limensi«)n:j:()t'  a  vo-ssel  (Jo),  turbulent 
cli^acteristics  pf  ^a  Hoty  j  =  1.2^...)..  etiu<litj()u.'»  of  iyiitiou  and  power  of  tlie  source  oi 
igmtibn:(^^-)i  Paraiiietei:  ids  of  the  tjtp.e  ()f  pdjaiueters  tliat.  lie  tiieasufed  raucli  easier  tliaii 
liCSiculatedi  Tliis^valuirwill  lnv.(Hie  of  tlie  parameters  of  our  niodel. 


3:i;-..fcc)adihg<:.bf;.u%alifeih:cti.sfe:blrleHa!iTati<nia)fca-iiuxtufe; 

Tlie  Regime  2^0.1  takes  phice  in  ctise  dimensions  of  a  vessel  . are  rtither  small: 


/■,)  <  i; 


(3:5) 


The  , velocities  ofliaine  . propagation  are  sui)spuic  aiid  rathei  small  in  this  c<use  compafed  with 
soiuiti velocity  fUid.it  ctm  b._e  assunied  that  presstire  luus  a  nearly  uniform  distfibutipii  inside  the 
vessel.  The^distribution  of  teinperature  will  not  be.  unifcnih  since  the  central  part  of  hiixtiife 
conrbusts  .underlow  pressure  conditious.and  is  coinpfessetl.theii, gradually  to  a  high  pressure  p„, , 
and-gassuear  the  walls,  is  coinpre.ssed  in.a(lvance  lietpre  the  coinbustibn  by  expanding  reaction 
.product^  of- internal  vtihunes.of  anixture;  The.  Kunperatlux!  rise  in- the  bust  case  will  lie  smaller. 
The  ^mpefature  ris.e  <lue  to  ajlitduitic  <;onipression;of  unlnunt  gas  can  be  deterniined  by  formula 


vn  -  > 


T  I  p  X— — 
-til 


(S.b) 


Po 


where  70  =  0;i/cU  -  the  ratio  of  specific  hesits. 

The  dependence  of  mass  rate  of.couibustioii  tipon  the  pressure  is  given  by  formula  (3.13  -  3.15] 


ih  =  pof  =  /Joei.i  — )  ■' , 

Po 


(3.7) 


where  71  is  the.  reaction  order.  The  formula  (3.7)  match(*s  most  of  experiments  when  n  =  1 
(3.14). 

Takitig  into  account  the,  compression  tif  the  unlnu  ned  mixture  in  front  of  the  flame  zone  and 
the  mass  burning  rate.  (3.7.)  it  is  po.ssilile  to  obtivin  formulae  [3.15]  describing  the  pressure  rise 
versus  time  in  case  of  flame  .sjuead  in  a  closed  v('.ssel: 

) .  for  spherictil  symmetry 


[Vivo 

/  P~~ 

\l 

/. 

Pm  -  1 

'(/P:  P  = 


Po 


(3.S) 


f  =  t  ■  iTr4l>oo„(P,„  -  1)/A/; 


2.  for  cvlindrical  symmetry 


7- 


(3.9) 


t  —  t  •  27r7‘()/ioi'n(P,„  —  l)/d7. 


The  solutions  of  these  equations  in  dimensionless  varitible.s  for  n  =  1.  7  =  1,4  are  shown 
in  Fig.3.1.  in  case  of  complete  comliustiou  of  a  mixture  the  upper  limit  of  the  integrals  (3.8). 
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{■3i9J;:is  P,,,.  .=  Ijf^  de.tennhiefl,  by  (d.S)-  '3.S).  T.3.'9),  nuike  k  pcwfjiidy  to  (letermiui 

tjn.=J.0,,,i},  Tiielipading  .(k',  thy  wall  i()i':r by  Regime  Xo.l  cau  !)y  dotenninyd  by  flonnuhi: 

■  ,  .  ...  ;/&!«  I  P^t-h  k-'/'^/,„-;  .  (3  10 

;  .'  '  ■  .  '  W  ;  I  fci  it  f  > 

inhere P:i:t)4nc\  Pi^  ,de.tenumed,;,fr(>iu,,.(3^S)i  (3’.9l.  (3.3). 

3.2L  rEbaflihgot’  walls  in  ca-soof  r<‘do*yion  of  liorinai  ciikonatioii.wavp 


,  ;r®^c,Bi6nxofku)i  luah  dbtptVtCU()jV  way^^^  place  in  case 


£  <  Co  —  / . 


(3.il) 


\^^here,/  is  thedeiigtU  of  unsteady  detonation  wave  piopagatitni  until  it  comes  to  a  .self-sustaining 
regime.  Tlieionset  ofcletonatipn  takes  plact*  at  c  =  t  and  the  pressure  iii  the  ye.ssel  at  the  time 
h|fld)ekl#ermmecld3y-fe^^^ 


Co  >■  Pill  —  I  • 


(3.121 


fpr.  the'case  tif.  .spherical  symmetry,. and 


^  =  ’l  -  P 

Cl) 


_iA.. -Pii 


(3.13) 


:f6r  tU6  ckse  of  cyliiidnt'fii  •syihiuettw. 

Corresponding  time moments  ti,  can  be  determined  from  (3.S),  (3.9).  The  time  of  reflection 
of  the  detonation  wave  frhm  the  walls  (*  is  estimated  liy  formula 


r  =  //.+ 


Co  **  L 


(3.14) 


where  D  is  the  detoiiation  velocity  for  the  mixture: 

V^2(7'  -  l)Q-r  \/2{‘r  -  l  )Q  +  4*;po/p»  , 


4t  <  =  t»  pre.ssure  on  the  wall  goes  ijp  to  /»,.  and  then  in  time  At  it  lowers  tlown  to  the  value 
PrD>  where  At  characterizes  the  reaction  time  in  the  detonation  wave  (3.1G): 


(3.16) 


where  E  -  activation  energy  of  the  chemical  r(>action.  P  -  universal  constant  for  gases.  p> ,  Ti  - 
pressure  and  temperature  Ijehind  a  leaiiing  shock  wav(>  prop«»giiting  with  tlie  velocity  D: 


D 


The  y^ue-oCAt  for  hvflrofteuroxYseu  iiiixturcs  is  ~  10“  ’  ^  m»l  for  flihite  niixturps  ueur 
fii|5defmiatipuXlimiU'i6-A’aiir^^^^^ 

iP;fessuresdh ‘thS-ueflecteilxiie^fomitw  yavp-aip  <Uf.tcruuji:!<l  by  fonluUaoi 
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(3.19) 


for  the-  leading  .dibck  wayft,- 
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(3;20) 


for  Uie;end-6ff  reaCiioivzm^^  ■ 

Tlie  unknown,  values  ()f  inessure  pos  d<*nsity  />p.  velocity  up-  at  the  endof  reaction  zone 
4)eiunclbthe;  detbnatiorii'M'aye  are-  determined  by  forimtliu*: 


up.  ^  ,  ,4.  „  <J^y 

<(iv  7  + 1  (hi  D 


(3.21) 

(3.22) 

(3.23) 

po  (tu  Pu 

The  pressure  near  the  wall  prp  decreases  exponentially  \inder  the  iiiHuence  of  rarefaction  wave 
,  that  follows  the  detonation  one.  The  characteristic  time  of  pre.sstue  decrease  is  that  necessary 
for  the  last  characteristicof  the  rarefaction  wave  to  rc'.ach  the  wall:. 


7  4^  1  </() 

IIP  ^  .  -D  , 

po  D  ~  i/p 

Pd  _  ^  ^  'iV.Pn.  ^  j  _j.  Po.fii'-D 

(<5 


r  = 


(3.24) 


where  Ar  -  the  distance  between  this  characteristic  an«l  the  detoiiation  wave.  «*,•  -  velocity  ot 
sound,  on  this  characteristic  (it^  =  0).  The  .solution  obtained  in  part  2  .shows  that  the  ratio 
^  =  itc  is  a  constaiit  value  verv  close  to  0;5 

ro 


Aj-  =  /.>/',)  s;  O.uco. 

The  value  of  (4.  can  be  determined  approximately  l)y  formula: 

.^+1  »ii3-7 


(Ik  »= 


D 
7  + 


1'  2  D-  2  ’’ 


(3.25) 


(3.2G) 


For  strong  detonation  waves  in  ga.se.s  \uider  normal  conditions  D  ~  2500  —  3000  m/s,  and 
flo  ~  300  m/s.  That  gives  the.  ratio  (illD~  ~  0.01.  and  formula  (3.26)  can  be  .simplified 


D 


(Ik 


Finally  characteristic  time  r  is  determined  as  follows 


3j3.^^.\hiaiiiic^i()a(iina-  of:tlH^Avaii  BBT-xouc. 

Tliis  type  of  ieHccubii  talws  piacf;  hi  . 


*‘i)  L.K  ;'i).  {3.2Sl 

The  transitioii.proc.ess  starts  after  r  =  £,au<l  thtit  is  why  siuiilai  to  Regime  K’o.2  for  the  present 
regime  ^e  shall  de.terinhie  tiiue  t\  pressure  from  the  solutions  (3.S).  (3.9).  Investigations  of 
the  pDT'process  (3.i>3.11.  Preliminary  Report  1]  sliuw  that  the  iiitensity  of  flelonutipn  wave 
increases  to  that  of  an  .overdriven  detonation  and  then  deeretisi-s  to  the  intensity  of  nonnai 
detonation.  That  is  why  tlu;  length. <)f  tin*  DDT-ztuu'  /  is  the  sum= of  nvti  lengths: 

/  =  .  (3i29) 

l\  T  the  length  of  forihathm,  of  aii  overdfiveii  dtitomitioin  h  -  the  huigth  of  slowing  down  the 
pverdfiveu  cletpiiatitin.to  {pself-siisttiinihg  nipde; 

The  ratio. of; these  leiigths  /f/j  ~  ()i()5  -r 0.5. 

To  begin  wifh.let  lis  fegiiard  the  .situation- 


t'o  /  <1  i'o  "  /i . 


{3.30i 


For  t  =  t*  pressure  on  the  wtdl  ri.ses  to  //,.  and  then  deert'tises  up  to  .p',^  during  the  time 
At',  The  pressvlre  beliind  the  detonation  wave  is  (h'U'rmiiusl  taking  into  aeeount  that  the  wave 
propagates  in  an  iiiitially  (•omprtissed  media 


(i  =  JliE  _ 

i>*  +  I  ai  7  +  r 


(3.31) 


where  p*  =  p(r);  «,  =  •  Detonation  velocity  D  di'termined  l'>y  formula  (3.15)  also 

depends  on  initial  paranu'ters  of  mixture,  hut  simple  estimations  show  rlnit  the  deviation  of 
detonation  velocity  D  from  tlie  initial  valu(>  is  negligibly  small. 

The  velocity  of  the  leading  shock  ot  .stu’h  tm  intensity  on  entering  the  unct 


will  be 


where 


unct  tnijnessed  mix t  un 


-  l>n 

.  6. 

(3.32) 

V  I'-i 

-  l>n 

I>n 

l>'>  _ 

£L  4.  i±JL 

/•'j  ^  -.-1 

I'n 

(3.33) 

p* 

-,-i  ii'.  ^  ^ 

The  Hugoniot  curve  gives  two  solutions  for  the  detomition  velocity  D'  >  D  determined  by 
(3.32).  These  solutions  correspond  to  strong  (overdriven)  ;md  weak  tletonation  modes  as  it  was 
discussed  in  Part  2.  Gas  velocity  and  pre.ssure  in  the  eml  of  reaction  zone  are  determinctl  liy 
formulae: 
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7+1 
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LD'  </„y>  2(7-^ -DQl 
run  D'’  u-  ]• 

(3.34) 

Z'l)  _  1  . 

—  It  ,  . 

(3.35) 

The  solution  with  sign  -f  in  (3.34)  corre.spon<ls  tt»  strong  detomition  and  it  must  be  used  for 
this  case. 
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'  \yaU,aij.<ir  tlM-  r(.*ti«'fti(>u  <jt  juioxiadrivcu  ('u*tpnati()u  wave  are  flet^nuiiicil 
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yi>§'4®PJ^6^se’pf,,i)rcsHUlv  tnnu ;//,,  up.  rak(isii>lac<'  (hiritisr  rinu^  niter.vai  At': 
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(l3S) 
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For  the  ca.se 
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(l40) 


(3i4h) 


pai’ameters  in  the  detonaticm  wave  and  at  the  end  of  reaction  zone  can  Ije  dp.scribed  by 
appfpxTmate  formiilae: 

/'^  =  -04m-/'m)(— ^  +  l)'+/4„, 
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P% - 'ip Dm  ~  Pm){  ■■,  --■  +  1)’  +  l>'p,„, 
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(3.42) 

(3.43) 

(3.44) 


Pirn  ~  /4 1 .l/pm  “ /’/j!i.=ro-J, : 

where  functions  //,  and  //q  are  dcitermihed  by  (3.31).  (3.35). 

The  wave,  velocity  can  be.  determined  by  (3.32).  (3.33).  the  velocity  of  -va-s  behind  the  wavt 
h'p  can  be  determined  a.s  ftillow.s 


iP'h  , , 
«o  iD"^p„ 


(3.43) 
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The  pre.ssure.s  on  the  wall  after  the  reflection  of  detomition  wave  behind  the  reflected  shock 
wave  jV;  and  at  the  end  of  reaction  zone  i>',[p  can  lie  iletermined  from  (3.30).  (3.37),  (3.3S).  Time  |' 
interval  At"  is  determined  from  (3.39),  (3.40). 

For  the  last  case  (r^  —  l\  <  L  <  co)  transition  proce.s.s  .starts  near  the  wall  and  is  not  finished 
by  the  time  of  reflection.  In  fact,  the  new  born  wtivt'  tloes  not  enter  fre.sh  mixture.  The  flow 
disturbances  during  the  stage  of  low  velociw'  flame  proptigiition  have  entmgli  time  to  reach  the 
wall.  That  is  the  reason  we  shtill  consider  the  wall  loading  for  the  times  f  <f*  to  be  the  same  ji; 
a.s  for  the  Regime  No.l. 


Ic 


j?ari  ,4,  Matijematical  niodeiiiug  of  spiiet'icni  .tuef  taiiks  fragmentation  under  tlie 
actipu  4pf  dyiian. 'c  "inteiaiai  :pi\essur^ 

To  analyze  caiisas  and  cousemuincch of  expipsioiis of  fncl  rankb  of  .laht  nx-kot  stages  oir  near- 
earth  orbits  there  are  first  experinunitson  fragnxcntaticin  of  thin-walled  shell  structures  under 
the  action  of! internal  explpsiph  (4;ij  iind  alscratUnnpts  to  solve  the  problem  thepreticiilly.  In 
[.4.2,-4-3).  the’  siinplest  matheihatical  inodids  are  suggested  rluit  all()w  to  calculate  ait  avenig«* 
number  of  fragments,  which  tire  the  result  of  iiretdmi)  of  thin  cyliiidrical  and  spherical  shells 
iu.exp.losipn,  and  their  initial  velocity  (if  scattering  provided  that  all  fragineiits  have  the  saiiie 
mass.  This  ptu't  deals  with  niiinencaj modelling  of  thin-wtdled  spheri<vil  slielLs  fragmehtatioii 
under  ihe  actipnad' dyntunic  internal  loading  <lescrii>e<l  in  Ptirt  3  taking  accoiiiit  of  fragments' 
distributioh.in  ief  nis  ofdhass; 

§4.1.  Problem  statement  dii  cleformihg  a  shelli 

The  follbwing  .assumptions  on  sliidl  gemiietry  and  a  charticter  of  its  deforming  and  breakup 
aremade: 

1.  The  shell, is  thin:  hjv  <  1  (,  h  -  the  thickne.ss.  r  -  the  sin'll  radius). 

2.  The  effect  of  ati  internal  dynamic  load  is  nuxh'lled  by  pressure  /t  =  i>(i)  depending  oil  timi? 
and -being  uniformly  distributed  along  the  internal  surface  of  tht'  shidl.  Charaideristic  time  of 
the. load  action  f  »  /i/ou  .{(in  is  the  .souiid  speed  hi  the  sladl). 

3.  The,  shell  material  is  considered  to  lie  ela.si(>viscoplastic  and  proci’ss  of  its  deforming  is 
adiabatic.' 

4.  As  a  condition  of  tla*  shell  breakup  lieginning  the  I'litropic  critiwion  of  a  limit  specific 
dissipation  (4.4)  is  assumed. 

5.  The  shell  breakup  is  assumed  to  occur  as  a  rt'sult  of  the  action  of  tensile  ring  stresses  at 
the  expense  of  consumption  of  elastic  eni*rgy  accumulated  in  the  shell  up  to  time  t  =  t,  of  the 
breakup  begiiming;  the  work  of  (‘xternal  forci's  in  lireakup  tiim*  is  neglected;  spalution  fractures 
aren’t  considered  (4.2.  4.3,  4.5). 

Due  to  the  first  three  assumptions  the  proldem  on  di'forming  a  thin  shell  can  be  considered 
as  the  one-dimensional  spherical  oiii'. 

Theii  the  momentum  e(piation  has  the  form; 

d>=  a!:!! -oil 
li  ~  r  ' 

where  p  •  demsity.  v  -  radial  velocity.  /■  -  currinit  value  of  the  .shell  radius,  a#  -  ring  stress  (the 
averaged  stress  over  the  shell  thickne.ss).  a  dot  over  a  .symbol  means  a  material  derivative  with 
respect  to  time. 

The  rate  of  ring  deformation  is  determined  l>y: 

a 

other  deformations  are  alisent  due  to  the  .sludl  thimie.ss. 

The  equation  of  mass  conservation  has  the  form  /)///  =  -'Hu  from  where  we  have; 

P  =  /»()exp{-2-//). 


where  pa  is  the  initial  density. 
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Typeset  lay 


,,  Stateequationsof  <4<K'jtpvis(!p|)iasti«\im]it,<'vi,ai  \vluch.  taking  arcotnii 
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;|o^&s^|ieriya^  cajf  . .  ;  ,  .  , 
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;l|ere  ./r  -;  shear,  modulus  //  -  dyiiamir  A’iscdsTty  of  inatcrijil.  ./u  -  >tatic  limit  of  elasticity  iii 
siinpy  tension V •  tlie  itiiit  <H,<jayysid<!  functipii.  HcrcuvifU  the  stn>.ss  tcnscu'  cr,  ,  is  divided 
dntq  sphericai.<T^),j;  andxleyiator  Sij^niftials  ami  it  is  assiunetl  that.iiites  c)t'  dftjforinatipns  can  lie 
xiiyidediihtp  '.tUe  rates  =of  elastic  and, plastic  ones  snul  a  piastic  How  isdncoihpressible: 

‘  i  ■  2^j'  ^'s{’i  =  '(). 

Specific  (l^er  niass  unit;)  elastic  tuiergy  E  and';inechan|cal  dissipafipn  D  caii  he  deterniined  by 
intei^ation  of' the -forhiiilmu 


*  —  T- - */») 


£)  =:  0  =Lrl'' 


.Rates  of  elastic  amfpiastic  defonuations  are  determinetl  i)y  the  fonnulaei 


it  “4  ijl  ii,  4.  —  41' 


2/<  3 

,  ColciUatiqu  of  the  shell  teniperature. in.  deftn-iuing  for  fidiahatic  approximation  can  be  per- 
fohhed  ii)y  the  folhiVving  pna'edure: 

lii'irf  +  2o,.5«r  =  0S«4''^. 

Herein  c„  is  the  heat  capjicitv  at  constant  stresses.  n„  is  the  cot*Hicient  of  volumetric  expansion 

(4.4 

Criterion  of  the  shell  breakup  Ijeginning  is  the  entropic  criterion  t)i  ti  limit  specific  dissipation 
that  for  the  medium  model  considered  is  redtu'cd  to  mechtmictil  dissipation: 


where  D,  i.s  the  constant  of  limit  specific  dissipation  determined  with  using  the  e.xperiments  on 
spalation  fracture  iii  a  phme  c<)llision  of  pl.-ites  (4.4). 

§4.2.  Calculation  of  fragments’  number. 

Fragments’  number  obtained  in  shell bretdaip  ctm  la>  found  frcau  the  balance  of  elastic  energy 
and  work  for  breaking  off  a  ipaterial,  To  describe  fragments'  distribution  in  terms  of  mass  in 
explosive  breakiitg  of  shells  the  Weilmll  distibution  is  most  oftcui  used  that  is  the  special  cti.se  of 
general  probability  distributions  (4.7.  4.S): 


.Y(<  ,»)  =  _,.xp(-(  — )■')). 

///, 

/8 


(4.1) 


Here  .iY(;g  -tu );  is  the  uiniiber  ()f  .fru'^uunits  with  iaass  ifss  rluiii  in:  .Vj  is  rhe  total  unuiijer 
ofifragiiients  with  mass  uioie  than  ()-(a  tWoretiral  <u)ijsr;mti;  /;/,  -  rhe cijiinictefistic  hiass  oi 
rtistnlntiihiiiuY?  the?iu.Ucx  <>f  frasinehtatitni.  uuality. 

;Fpr.rA.>,  i  tlistrUniti()u444)  is  vuh  with  the  iiKnle  -V/  —  for  A  =  1  this 

-ciistributiqnibecpmesithe  expbueutiai  .()ue:,tbr  ()  <  A  <  1  it  has  au  asjhapcote  ireiag  the  ohbiiate 
■0:is.  yalue:gfc.:Vo,is-fihite44  A-  >-,L 

As,  shpwn».in  ;(.4J.  .4.Si,  the.  luiimpchii  (.listributuHi,  (4.1)  4, escribes  satisfactorily  spectra  of 
breakup  pt metallic  .cyHaclfical  sheils.  of,  a  .aw'tiimii-.sizcil-  tluckness  (li /?■  ~  0.1 ).  Better  results 
cahvbesgiven  ;by  the  :tworaiotlul;liype.rweiljul  <listril.>utioa  sugfte.sted  ia  J4.71  which  iaclutles  two 
morphological  cpUectipus:  large  fragiueats  cpataiaiag  Ixah  shijll  stifftices  (  exfcraiil  aiicl  iuteriial) 
ahHiaccpmpahying.-.sinallei^fragmei  cpataiiaug  .one  of  tlie  stirface.s. 

For  thexase  c.pn.sidefefl  siielis.are  very  tliia;  h/.r  %  O.OOl  (4.1).  Therefore,  evidently,  \ye  caii 
restrict  ourselves  to  luamotlal  .distributipa  ^4.1)  aad  a.ssmni*  that  till  fragments  contain  both  an 
inteni^  surface  and, an  external  surface,  of  the  shell. that  is  <'oiifinhed  iuclirt’cily  by  experiments 
(4ii)  conducted  for  suebthiinshells. 

Letran.Mea  of  an  initial  externaLLagrangian  surface  of  a  fraguu;nt  be  eijiial  to  .s,  its  interiial 
.surface  be.^.so  equal  to.  .s  (due  to  the  .shell thinne.ss),  and  aretrof  its  lateral  surface  -  2ph,  where  p 
is  the  semipefimeter  of  the  contour  s.  The  fragment  iha.ss  is  m  =  /)(,//. .s.  therehire  distfibutipn  in 
terms  of  ma.'js  (4.1 )  can  lie  presented  in  the  tbrm  of  di.striliutiou  in  terms  of  areas  .s  of  frcigments: 

A'(<  .s)  =  A'„[l  -  ,.xp(-(-l)‘')|.  (4i2) 

•> 

where  s,  =  ru./poh  h  the  characteristic  area  (if  the  friigment. 

Number  {»f  fragments  with  aims  s'  <  .n  <  s"  .will  be  equal  to: 

j\’(.s'  <  .s  <  a")  =  At,, j  ^ 

Suppose  now  that  from  the.  tcital  spectrum  of  the  shell  fragments  with  aretis  0  <  .s  <  Ttds  (do  is 
the  initial  diameter  of  the  spherical  shell)  one  can  (extract  K  groups  of  fragments  with  area's  .sq, 
■S2, Sk'  Siniu  <  ■si  <  >■’2  <  .<  •'''/v'  <  qini/aud  s,„iir  i'le  certain  minimal  and  maximal 

areas  of  fragments.  Let  all  fragments  with  areas  .s„,,„  <  .s  <  0.5(.‘ii  -f  )  prove  to  be  witliin 
the  group  of  fragments  with  urea  .S|.  all  fragments  with  art'as  ().o(.'‘i  -r  .'•<)  <  s  <  0. 5(^-2  +  'i) 
prove  to  be  within  the  »r()up  with  ari'a  .nj.  etc.  uj)  to  the  grotip  with  .uea  wdierein  there  are 
fragments  with  iireas  O.o(sr— i  -i-sf.)  <  .^  < 

Further  in  place  of  distribution  (4.2)  the  following  distribution  will  lx*  used: 


iY(<  =  A’o[l  -  exp(-{ 


•'*  “*  •''min 


)  )].  •'‘iiiiii  < 


ft 


Then  the  number  of  fragments  of  the  group  sj  {j  =  1.2 . K)  will  he  etpial  to 

, (lo(.’'^_i  4*  Sj )  “  •' iMiii  \ 1 


A^i  =  A'o  ( i^j  , ) ,  ij  =  exp  ( -  ( ■ 


■r). 


(4.3) 


(4.4) 


j  =  1:2.  ....A'. 

Here  so  =  2s, ,,,-,1  ~  ■'*i.  .^/c+i  =  2s,„„j.  - 

System  (4.4)  of  K  equations  for  cjdculation  of  the  fnigmt*uts'  nuniljer  of  the  groups  can  be 
complemented  by  the  following  two  ecjuations; 

A" 


2=1 


»  £ 


fl 


}  . 


•  'iL 


h 

ft 
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.w>n.«  4*..- 


l<  . 


(4iG.)= 


-3=,1 


;Eqiiati()i\  (4..o,)/ni(fms;tlua  s\imiuurizc<l,ar<'a  of  (•xrrrii;il  .siufar<‘.s  <)f  tlio  ti-ai^uients  is  exactly 
equal  to  tiie  41*^1!  afeju  iuul  «;(iuatiou,(4.G)  uaiaus  that  tlie  eliistic  eii(‘rgy  r</,',/i/A)E.  fEi  ~  I 

accuuiulaterl  ill  till*  slitiil  .by  tiine  f.  ==  f,  si).<*n<ls  fof  formation  (if  Ine'akiip  suiitaies.  The  coeificieiit 
I'E  ^  t)  oficlaLstie  (".;'tu;y  cohsimiiitiou  forifonuatiou  of  I )n*akui)  suifac(*!j  is  iu'tnxiuccG 
into,  equation  (4.6):  Furtlit*!'  \vas  {\.ssuin(*<l  in  .calculatitius  that  all  tdastic  eiangy  spends  lot 

fgnuatipn.of  bietdiup  (ial.  hr;  =  lb  .Hennn  7  is  the  speciiic  Jeneigy.  coiisunied  foi 

fpniiatiompfdlie  fret*  surface  unit,.  //,  is  the  .senliiierinlitter  (if  a  fragineiit  of  the^/Vih  gioiqi. 

Gil  u.sin|  (4;5,);  thetfollowing  eqtiatitni  can  be  (ieiivethfrom  (4.6): 


if 


' '1-  t 

]=l 


(4.7) 


Tq.find.(27v  +  1)  uulnipwn  variabl(*.s,p^i  S’j.  .Vo  tluTe  are  oiily  (A  +2)  (Hpuitions  (4.4),  (4.5) 
ahd:(4;7)(  S(*lec.t  thti  ftilhmduii;  particiihir  stihituiiis  of  eqiiiitioU  (4.7): 


V'li  At /.'/■;  7 

l>,  =  — : - -s,.  1  =  1.2. 


.K. 


(4.S) 


Phy.sical  nieaniug  of  (^iS)  consists  in  the  followiii";.  One-half  of  tnieiAy  necessary  tor  forination 
of  breakup  surfaces  htouiuI  each  ./-th  frasnieut  with  area  .Sj  and  mass  iiij  is  extracted  from  mass 
mj  contained  insidt*  tlu?  frasimnit.  and  the  .secoiid  half  -  from  the  outside  (froiii  nei,c;hborin.ui 
fra,9iieiit.s).  This  is  jqiparently  a  reti.soualili'  assumption. 

The  iust  of  (K  -H)  ('(puitious  (4.4).  (4.G)  eau  l»e  easily  .solv(*d  now.  ,\t  first  we  find: 


.Vo  = 


-t/ii 


S!=i  •'*i(-A  ”  h+i ) 


(4.9) 


and  then  according  so  (1.4)  all  Xj  can  be  determined. 

Notice  that  if  wejiitroduce  the  non-(lini(*usional  coefficient  of  shapi'  h'j  =  'jlp'j  then  we 

will  obtain:  _  , 

/,  =  [M^EEI-,  (4.10) 


For  plain*  figures  the  shape  coefficient  changes  within  tla-  limitj 

U</.-,  <1/-. 


(4.11) 


provided  that  its  maximal  value  is  the  one  for  ti  circle.  .4s  iollows  from  (4.10),  fragments  of  a 
greater  area  me  le.s.s compact,  that  is  verified  !)y  experina'iits  on  cylinders  (4.S).  Small  fragments 
can  be  considered  as  plane  ones,  therefore  the  following  ivstriction  from  below  on  the  fragments 
area  is  derived  from  (4.11): 

ii»t,hr: 

that  is  necessary  to  take  into  act’ount  when  is  Ix-ing  cho.sen. 

It  should  be  noteitl  that  if  the  characteristic  liuetir  dimension  of  a  fragnumt  Ij  =  2.<j/pj  i.-' 
introduced  then  it  will  prove  to  be  tlmt 


puE^hi; 
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1,65  the. chi^actfnis.tic.;li  climeu.'^ion  the  hame.aiKl  coincifleut  with 

oUie  pne-^^pbt^iued  at  solution  oi  the  proiilein  when  all  t’ia“UKnits  are  assuiiied  to  be  the  sauie 
tsplutidh  )■.- 

|4s3i  The  guiciUig  paraiWeters 

■Fpr  the  ctHisidered  probleiii  pu, shell  fragnieututipu  there  is  a  iimuber  of  ptiranieters  which  caii 
.be^d,iyicled  into  two  parts:  paiaiueter-  idiaracrenziug  ilie  shell  material  (pu,  //.  ;/,  (>.  o,.,  Ju.  7. 
,l)i)vaudi?aramet6rs  tm  which  a  ciistributicia  <)f  hagiuents  in  tefhis  of  nni.ss  l»u(vS«).  A)  dopphds. 

For  the  kiiowii  poiistruetional  iiuiierials  such  jiaraiueters  as  density  under  the  normal  con¬ 
ditions  ppi  shear  inoduhis  p.  heat  (capacity  at  constant  stresses  (>.  coelficieiit  of  vohuhetfic 
exR^eion  O',,,  static  elasticity  modulus  in  simple  teiisiou  ./„  are  table  vtilues;  It  is  naturally 
that  in  deforming  process  with  changing,  tc'inperatun^  T.  .(huisity  p.  pressure,  growth  of  plastic 
deformations  these  pafaineters  also  change;  Flow  limi,t  ./,)  and  shear  modulus  /i  are  assumed  to 
change. iu;  tin}  defprining, process  according  to  the  Steinljerg-Guiium  model  (4.9): 


=  .JUil  4-  -'A(t  -  To)). 


p  =  p„(i-/w(^)'/‘-')((r-r„)). 

I' 

Here  Jqi  /<o  are  values  of  parameters  under  tin*  normal  conditions.  s',\  =  2|i{J{  is  the  intensity  of 
plastic  deformations:  /h  n,  b.  J,„,ijr  are  constants  of  tin*  matefiid.  These  constants  for  many 
mateiials  are  pre,sented  in  (4.9). 

To  determine  dynamic  viscosity  of  matefial  >;  is  more  complicated.  However,  tins  parameter 
is  also  presented  for  many  materials  in  the  papers  using  the  motlel  of  a  .solid  rleformable  body 
quick-response  to  deformation  rate  (f.e.  (4.0)).. 

Paraineter  7  (energy  necessary  for  formation  of  tin*  brt'akup  surfai’ci  unit)  is  also  shown  in  a 
number  of  papers  (see,  for  example,  (4.7,  4,8)). 

Now  consider  briehy  determination  of  parameter  D,  being  limiting  speciKc  flissipation  (crite¬ 
rion  of  the  Ineakup  beginning).  .A.-,  it  was  noted  in  !i4.1.  parameter  D,  can  lie  determined  using 
the  experiments  on  spalation  lueakuj)  in  plam*  collision  of  plates  anti  numerical  modelling  of 
this  process.  It  can  la*  done  by  the  followiiig  way. 

In  experiments  a  plane  collision  of  two  thin  phites  occurs;  an  imptictor  :ind  a  target  (see,  for 
example,  typical  experiments  (4.10)'..  ..^s  a  result  of  shock  interaction  of  compression  and  tension 
waves,  in  the  target  under  the  certain  conditions  the  region  of  tensile  stresses  expressed  brieliy 
appears,  these  stresses  lead  finally  to  formation  of  a  spahilion  cnick  and  in  a  number  of  cases  also 
lead  to  a  complete  breakup  of  the.  target  with  withdrawal  of  ti  so-called  spahition  plate.  Spalation 
fonnation  becomes  ajjparent  in  the  grtiph  of  depemlenci*  of  the  velocity  of  the  btick  target  surface 
on  time  w  —  t  registered  in  the  experiment.  Such  typical  experimental  dependence  xo-t  is  shown 
in  Figure  4.1  for  the  ca.se  of  collision  of  tdumiuium  imptictor  with  titanium  ttirget  for  two  collision 
velocities  Vo  =  GGOm/s  and  Vo  =  1900m/s;  stdid  curve  corre.sponds  to  experiment  (4.10),  dashed 
curve  -  to  computed  results  [4.4].  In  a  ptu-alh*!  way  numerictil  ;.i()delling  of  collision  of  phites 
is  conducted  using  computations  with  explicit  separatiouof  surfaces  of  spalation  breakup  of  the 
target  ( in  detail  see  (4.4)).  In  one  of  the  exp«*riment.s  by  means  of  vtuying  the  model  parameters 
(  in  particular,  the  parametr  D.)  the  be.st  coincidence  of  computed  results  with  experiment  is 
acliieved,  and  the  rest  of  experiments  are  used  for  control  of  parameter  selection  quality.  Thus, 
for  titanium  in  (4.4)  the  value  Z),  =  75  k.J/kg  is  chosen  th.it  tillows  to  descrilie  (juite  satisfactorily 
a  number  of  experiments. 


Gh^acter  of -fe  (ti-’''rir)uru>iL;iii. ftTigs  i>y  two  j)aniiuc*t(n's; 

ci^acteristic  an*u  .s,.  ,iiu<l  iu/Urator  of  ci'usluus,<iu:ility  slioliUl  ‘ 

on.the'basis.cif  tho  experiuu'utaVdata.  XliVis.  lor  rliis  oljjcfrjVc  ('xpyiiiiii’iitS  (4,1]  \yii|  bo  tisofU 


5.  j^^meiits  acceleratidii  by  gas  cidud  propagating  into,  vacuum. 

iiii?  part  of  oiu'  investigations  is  denoted  to-obtaining  the  rate  of  final  velocity  of  fragment^ 
fonpeci  by  explosion  of  the  fuel,  tank  alter  breaking  up  its  walk  U>iug  soiiie  assuiuptious,  we 
reduce  the.  lutniber  of  dimensiimless  goviuliing  j)araui(?ters  of  the  jnoblem  to  4:  gah 
adiabatic  con.stant,  relation  of  gas  and  fragiueut  s  <l(.'usities.  initial  vtdocity  of  a  fragment  relati'd 
to  the  Rayleigh  sound  velocity  and  the  g<*u«*rali2e<l  drag  (•o('ffi(’i('nt.  \Ve  made  computations  for 
sphefical.and  cylindrical  cases  investigating  tire  vast  diapasone  of  jiiirameters.  obtained  arrays 
of  the  filial  fragments  velocities  and  then  \ve  managed  to  ol)tain  high  precision  approximation 
fonnvtla  connecting'  the  fintil  dimeusioiihtss  vi'locity  of  a  fragment  with  the  4  pai'ameters  liieii- 
tioiied  above.  To  obtain  then  the  real  velocity  oue  lias  only  to  inultijily  it  by  the  Rayleigh  sound 
velocity  ctilciiiated  for  the.  ihitiaf.ptu-ameters  in  the.  gas  cloud, 

5.1.  Probleni  statement  aiid  assumptions. 

We  assvuue  that,  the  initial  gas  cloud  is  uniform  ami  has  a  .synmu'.trical  shape  (spherical 
or  cylindrical).  For  it  propngtites  into  vacuum  tin*  disturbances  in  it  decretuse  in  the  process 
and  usually  after  refiectihg  from  the*  wall  propagate  to  the  center  not  tirtecting  the  fragments 
originated  from  the  liroken  wall  so  w(*  assume  that  the  uouuuiformities  of  jne.ssure  in  tiie  inititil 
ga.s  can  not  affect  heavily  on  the  rate  of  the  Hnal  results.  Howevi'r  we  must  jit  first  olitaiu  the 
average  rate  of  pressure  ttiid  temperature  in  tlu'  initial  gas.  We  do  it  according  to  the  following 
fonnulae: 

«  li  H  n 

P^=  J  r^vilri^j  =  J  r'Tdr(  /  (o.l) 

e  0  (I 

where  u  equals  to  1  for  cylindrical  case  and  2  for  sphtnieal  one. 

We  assume  the  initial  velocities  in  the  gas  eiiual  to  zero. 

The  main  assumptions  concerning  the  fragments  are  the  following:  we  consider  tluvt  the 
fragments  affect  on  th<^  gas  much  less  than  the  gys  affects  ou  tlu'  fragments  .so  we  neglect  tin* 
first  effect.  Then  we  neglect  the  fragments  infiueuce  on  each  other.  So  we  assume  the  fragmeuts 
independent  and  not  affecting  the  gius  motion. 

With  the.se  assumptions  we  obtain  the  followin.g  mathematical  statement  of  the  problem. 

We  will  use  the  Lagrange  s  tipjiroach  to  the  proldeiu  considering  the  imiss  of  gtis  from  the  c(Uitei 
or  the  axis  to  the  fiuid  locus  as  the  Liigrniigiau  coordinate’.  This  iipprotu  h  to  gas  dynamics  is 
widely  described  in  literature  (l.e*.  iu  (5.1)).  The  Ltigriuigiiiii  um.ss  coordimite  here  is  determined 
as: 


~  j  (5.2) 

II 

A.t  first  we  introduce  the  diineusiouless  coordinates  and  varitibles  Jiccortling  to  the  fonnulae: 

P  ~~  ^  TflT,  E  —  EnE,  fi  =  /tti/l)  /'  —  /?/’,  ~  5i)s.  t  ~  fuf,  c  —  (e.S) 

where  p,  T,  E,  p,  r.  .s,  /,  c  correspoud  to  pre.ssure.  tempeiiiture.  internal  energy,  diuisity. 
radius  (distance  from  the  axis  or  the  center).  Lagrangitm  mass  coordinate,  time  and  velocity  of 
the  gas.  The  scaling  rates  P(,  tiiul  To  are  considered  to  1h‘  the  average  pres.sure  and  temperature 
in  the  initial  gtis  correspondingly  determined  in  (5.1)  and  the  other  scaling  rates  are  determined 
by  the  sequence  of  formulae: 


Po 


Ei/Pg 
U  To  ■ 


r„ 


i;.’ 


(5.4) 
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Typeset  by  .4 


•EiiTejE;A  =  SSI  r,h«  nuivi-i>;il  i^«s  <«)iistaur.  !l  >  t.iu'  avcnt“(‘  molar  mass 

aiKlfi?o-is  tile  iiiitial  radiiis  of  tia-  !>;as  cloud. ,X)iie  can  >.  «•  from  (5.4)  that  the  velocity  scalimt 

rate  IS  th6  Rayleisli  stiiiiKl  vidocity  cakulat(*(ffAf  ihuiariias  iiaraiiietc-rs. 

Below  tye -shall  remove  tilde  upon  lettofs  in^thc  foiiimiae  a.ssiiminu;  all  the  parameters  hi  the 

g^  todi'e  tlliueiisioidess -if  not  specially  iiieutiiiiU'tl.  ^ 

We  iiaye  the  iiaiaiices  oftiutuss:  iiupiilse.aiul  euerf'y  in  t  he  .i^as  to»;ether  with  the  state  etpiatiou. 
theiutenittl  euersy-eepiatioii  aud  tiie  Euler,  s  ratlins  etpiaiioii  iii  the  L(e4i<tnu.iau  <  ooidiiiatt  ■  y  • 

'  '  '  Ot  '  Ok'  Ot  -Ot^r 


-  T  r  -  ,<  (5)  , 

,,  =  ,.r.  E=—.  . 

Here  7  is  the  ga.s  adiitbatic  mtio  lussumed  to  he  coasrtmt. 

the  Ijoxiudary  tumditioiis  are  determiued  iu  the  ceiitt'r  or  011  the  a.^is  (.<  =0)  and  tit  the  edge 
of  the  gas  cioud  (.s  ai)-  As  we  forset  that  the  gas  propaittites  into  viicunm  the  value  iit 
is  constant  iu>the  process.  The  value  of  determinethticcordiiigly  to  the  deternimatum  ol 

the -Lagrahgian  mass  variabh'  (5.2): 

1 

—  f  r"  mil- 
() 

So  we  have  the  followinii  bouiidavy  conditions  to  de.scrilie  gtis  propagiuion  into  vacuum; 

,S  =  ():  esO:  =  /' =  0- 

The  initial  conditions  tue  simple  itccortling  to  our  ii.ssumptions  and  the  way  of  determining 
flimensionless  values  (5.3): 

/=:0:  e  =  ().  p=l.  r  =  l. 

We  can  see  that  the  •solution  of  the  problem  in  -..eneriil  iu  dimensionless  lorili  depends  onl\ 
on  the  value  of  the  {idiahatic  coetticieut  *. 

A  fragment  in  the  gas  How  has  the  folltiwing  laws  of  motion  tin  dinientional  vtirialiles): 


where: 


S  =  Ecoso  +  (^  +  )£.-■/ tio. 

•iifi) 


C,  =  C,(o./).  f  =  L'/Z. 


Here  m  is  the  fragmeut'  nia-ss,  »  is  the  fragment’s  velocity.  E  is  the  area  of  one  side  of  die 
fragment,  6  is  the  frasaient',s  thickne.ss.  C,/  is  tin*  dia»  coefficient  (assuniing  the  velocities  high 
we  neglect  its  dependimce  on  the  Reynolds  uumlier).  S  is  the  eHective  area  facing  the  How.  o  is 
the  angle  of  fragment's iuientation  (zero  if  it  is  faced  perpendicular  to  the  How),  L  is  the  eH’ective 
.size  of  the  fragment.  /  is  the  slnipe  coefficient.  Here  we  ii.ssunie  the  Iragtnents  originating  tioni 
the  thm  membrane  will  of  the  fuel  tank  with  sizes  much  less  than  the  diameter  of  the  tiuik  and 
the  fragment’s  surface  «‘ur vat ure  radius  equal  to  /?(,. 


i 
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4 


4 
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Tlie  fragiiujiic  iii  the  How  ror.ntes  oi  course,  hut  we  .'.hail  a.s.siime  rhat  iu  its  rotation  the  eHective 
orientatiou  to  the.  How  o  keeps  its  value.  The  shape  coeHicieut  /  never  changes  in  the  process 
for  one  fragntent;  .so  we  obtain  that  tiie  value-ot  the  <*.'Cpre.ssion  {C,iS)/2  in  the  second  term  ot 
tlie;impulsh  equation  (5,.91.is  constaiit, 

The  initial  conditions  for  the  fragment  iti  dimensional  form,  a  re  the  following; 


if  =  0  :  11  =  (/„.  .(■  =  /?(,. 


(5.10) 


iisiiig  dimensionless  vtuiables  and  con.sidering  the  scaling  nit<‘  for  .c  ('qiial  to  i?o  and  the 
seeing  fate  f(jr  u  equal  to  Vd.  we  can  rewrite  the  laws  of.motion  (5.9)  aiid  the  initial  conditions 
(oTO.)  for  fragments  ( reiuoving  tilde ): 


It  ~  ~  “  ''H"  ” 


(h 

ili 


=  It. 


f  =  0: 


(/  =  V.  X  =  1. 


(5.11)- 


(5.12) 


The  dimensionless  numbers  0.  D  and  L  correspond  to  densities  ratio,  generalizetl  drag  coef- 
Hcient  and  dimensionless  initial  veloc.ity.  and  are  iletermined  accorcling  to  the  following: 


Q  =  —  Hi  Q  —  ^ 

>11  !>{  2E(') 


U  = 


"(I, 


(5.13) 


The  result  we  need  to  oljtain  is  the  value  of  fragment  's  velocity  tit  the.  infinite  time  from  the 
begining  of  the  process,  so  it  is  indepemlent  of  the  tinu'  scale.  We  will  note  it  Uf^,  (dimensionless) 
or  Uoo  (dimensional).  To  obtain  the  dimensional  vtdocity  wi*  must  multiply  tlimensionless  one 
by  Vo  (Rayleigh  sound  veloc.ity  for  the  itutiul  gas  ptinuneti-rs). 

So  using  our  assumptions  we  reduced  the  set  of  jtarsimc.iers  tifiecting  the  dimensionless  final 
result  to  7,  G,  B,  and  U .  Below  we  will  de.s<nibe  the  comptuational  model  obtaining  ai'rays  of 
Uoo  wid  finally  the  approximation  formula  taking  into  ac.cotmt  these  four  ptirameters. 

5.2t  Numerical  model  for  gas  cloud  propagation  and  fragments  motion. 

To  calculate  the  gas  motion  the  numerical  monlel  was  used  (Icscrilicd  in  (5.1).  We  used  the 

differential  network  ot  .V  =  513  luxles  /.•  =  i) . .Y  -  1.  ca.-fi  node  corresponding  to  a 

definite  value  of  Lagrangian  nuiss  <'oordinate.  The  network  was  ljuilt  uniform  relatively  to  the 
initial  state  of  Eulerian  radius  (non-uniform  ndatively  to  thi'  Lagrangian  coordinate).  We  used 
implicit  scheme  to  approximate  the  gas  dynamics  c(iuaiioji  and  a  "plain  iteration”  techniqvie  to 
solve  the  prolilem  cm  each  st<>p  of  time. 

Actually,  the  transfer  from  the  Ji-th  to  n  1-th  time  Itiyer  is  ol)tainecl  from  the  equations 
(5.5),  boundary  conditions  (5.7)  and  initial  c-onditions  (5.S)  and  looks  like: 
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•rf*' -iC, )  =  i).  '.-  =  1 . -v-i 


^-0.5(  +ei!)  =  ().  /.=() . .V-1 


'■r  ’  -  ' 


ypK  .p  ^  ^  ^  Q . Y  -  2 

-Et+Q.5(v'r  +K!)(l//'r'  -1M!)  =  0.  /.•  =  0 . .Y-2 
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t,  I  U'l-  '  T'.a-  -r(t^.)  )/-J- 

Hi-  = 


1/  ==  2. 


/.='() . N^l 


(/•r' -rnj:)/2.  //==!.'  , 

J«+ 1  ^  £..+l  ( .  ^  1  ^,»+ 1  ^  j»-i- 1  ^,»-r  I  y . ;Y  _  2 


(5.14) 


Here  the  value.s  of  iy,  /'t.  Ek  ‘iro  ("ilrulaicd  in  rhc  /.-rli  nodi’  .^A  of  the  differential  .sehenie. 
the  other  values  are  ealcuh.ted  between  .•»*,  and  ••'t-t-  !n-  =  -t  +  i  ~  •'t-  ”«  timesti^p  value. 

in  order  to  supress  oscilu.uon.s  we  use  aiuticial  '.i  umann  s  viscosity,  addiiif;,  it  to  pre.s.suie  in 
the  iihpvtise  e<[uation:  it  is  eaieulated  l»y  the  .seijuenee  oi  iminulae  between  ••'a  and  .•'k+\  aeconlina, 
to  [5il): 


,L  = 


e,. 


ii+i 


A-  . 


In. 


,,^.1  ^  5(7 --Ilf  '/.•<«• 

-Y-'  •  1  0.  il„  >  0. 


The  bouiKltuy  conditions  in  the  ditferentitd  scheme  an': 

i<  +  l  .  (i+l  _ n  —  o 

'()  - /'V-i  ~  "V-l  ~ 

The  initial  conditions:^ 

vl  =  0.  ijI  =1.  Tj!  =  1.  Et  =  1/(7  -  1 ).  =  0 . y  -  1 


(5.15) 


(5.1G) 


(5.17) 


The  algorithm  realising  tlie  scheme  (5.14-5.17)  is  descriln’d  in  (5.1).  We  u.sed  the  simplest 
techniqvie.s  of ’’plain  iteration":  it  retiuires  tht'  modifii'd  Courant  criterium  for  the  timestep  rate. 
Actually  we  used  the  criteritnn  taking  the  Neumann's  viscosity  and  possiliility  of  In'eaches  into 

account: 


k>{):  ().2(r;t|-i 

/.•  =  (): 


)/''A. 


where 


rtfc  =S  < 


I,  ,  .  ..  ,  ,  ,  1  »+i  .M+ii  ^  t 

Y  -f  I )  +  1  -  1)  ^  - ^ 

V  ^ /  +  u  +  /  -  1 )  i '  a  ''a+i)' 


(5.1S) 

To  calculate  the  fragments  motion  we  itseii  tiu’  "predict  or-coiiector  technitiues.  or  the  modified 
Euler  tedmiiiues  of  the  secontl  order  tipproxiimition  {5.2).  We  calculat(*d  the  fragments  motion 
on  eacli  step  of  time  after  ciilcuiating  gas  motion,  so  we  ttsed  the  same  rtite  of  time.step.  .A..s  the 
fragments  do  not  affect  gas  tind  etich  otherHn  our  model,  we  treated  <i  l)ig  numljer  of  fragments 
with  different  parameters  simulttineotisly. 

Actually  to  olitain  iuriivs  of  ti-o  we  examined  the  following  ditiptrsones  of  parameters:  G  = 
10-4  ^  10~^.  i?  =  10  -i-  1000.  C*  =  0  -T-  1.  Then*  were  7  viilues  of  parameters  taken  from  each 
diapasone  including  the  lowest  and  the  highest  value,  so  343  model  fragments  were  examined  in 
eacli  calculaticn.  .A.ctually  the  set  of  values  used  in  calcuhitions  looks  like: 

G,  =  {1'10-*.  2'10'-‘.'o'10-‘.  M0-‘.  2-10-*.  5'10-‘.  110"M- 
Bj  =  {10,  20.  50,  100.  200.  500.  1000). 
f;A=  (0,  1/G.  2/G.  3/G.  4/G.  5/G.  1). 
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The  calculatipn4ts(4rde4^‘"4tKM^-'4^‘>  rlir-valiic  of*  auil  pii;  tlu-  rypi'  (jt’ ii,e()nietry  (spiicncai 
or  .cylm4ricai)r  ‘’U  ''-.  Sc)  yv  laade  calailatipuh  for  =  {1.20;  1.23.  1.30.  1.40}  aiai  for 
botU^sgeometfies  <d)tmuiipi!;!ih-;(N  (-a.sc*  array  .of  ii^  valiuvs.  , 


K3i 


Such  a  |Mv  iuimliei'  (,)f  dinuuisiouless  i)aiiuucter.s  mako.s  it  possihh*  to  try  to  obtain  approxi¬ 
mation  iformviia  for  Urx,:  pf  tliis  ipfinuia  t-air.savf  calculations  rinu*.  This  fcnmuhi  needs  to 

i)e  precise ;eiumghdn  the  diapasone  of  parame,ters  .mehtioneU  above. 

To  obtain  the  formula  we  used  inininuil  sciuares  rcs-hnhiues..  however  slightly  incKlilied.  Let  ns 
ihy6rye4he,;folio\viiig'n()tjiti()ii; 


-y'lUi-  t 
"  bo  —  ' 


uT''  =  \V(,^.-,„„G.D.U). 


where  .the;first  term  is  the  calculated  array,  iuflextfs  I,  correspond  to  the  values  of  G,\Bj.Ui, 
uswi  ih  calculations,  and  the  second  oiie  is  the  approximative  function  we  need  to  olitnin.  The 
nhihimal  squares  techniciues  retpiires  us  to  minimize'  the  functional: 


0  It  It 


laU  /=U  t'=(| 


vj;^ 


.7.A- 


For  we  have  ctilculated  the  tirruy  V'  for  certain  vtilues  of  t  and  //.  then  wt*  neetl  to  suggest  ti 
set  of  Vr  formulae  for  these  values. 

We  suggest  the  following  formula  con.sistirtg  of  10  coelticients: 


C»GD  +  CiG  1 

C.iGD  +  CUG  +  Cr,L’  '  C,U 

1  -f  c\iGD  +  c\  G  -far  ^  1  -r  Cy  u 


We  have  solved  the  problem  of  <!•  value  minimization  for  the  values  of  7  iind  u  mentioned 
above.  As  the  mininmm  is  close  to  zero,  then  to  find  the  values  of  C„  array  with  maxinnil 
precisity  we.  actutilly  minimized  not  ‘I*  itself  luit: 

,  f  log»I>.  <I»  <  1. 

vp  =  { 

1<I*-1.  ‘I>  >  1. 


We  applied  the  heavy  Ijall  technup;es  to  obtain  the  minimum  of  'P  [3.2j.  The  techniques  ot 
gradient  .slide  [5.2]  here  is  less  aviiilable.  foi  the  l()-<limeusions  'P  hypersurftice  has  narrow  gaps 
and  therefore  the  gradient  .slide  retptires  much  more  steps  to  find  the  minimum.  In  order  to 
optimize  the  process  of  minimization  the  friction  coetticii'iit  was  different  for  different  flirections 
and  it  raised  automatically  with  (\scillatious  of  the  "baU". 

Our  results  are  shown  in  the  ttddes  1-2. 
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Tubl.H.  C(iiiffi(:u:nix  for  cjilinilritui I  tui.'tt:  (u  —  1  j. 

It  can  bo  jiotioed  mnu  the  fiil)l(*.s  1-2.  that  tlu'  prooisUy  i»t'  appioxiuiiitiou  i.s  lu,u,h  .so  that  tho 
rate  of  relativo  prooisjty  of  1%  in  tho  oxamiuotl  (liii[)aM)in‘  of  piiiaiuotcis  i>  i!,narauto(l  (actually 
the  precisity  is  much  hiu;hor). 

To  obtain  tho  approxinnition  foimula  takiuu,  *.  into  accouui  wo  use  tho  Lai!,raas!;iau  polynom. 
•SO  that  if  IT,„  i.s  calcc.lattnl  for  'i,,,  accor<liu»  to  tho  talilo  1  or  2.  tho  value  of  IT  for  7  i.s  calculated 
(see  (5.2)): 
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Tliis  approximation  formula  is  the  main  result  of  this  pait  of  out  report.  It  can  sig,nificantoly 
simplify  the  cidc.ulations.  To  obtain  tla-  dimousional  value  of  tho  fra!!,mout  "s  Bmil  velocity  one  ctui 
multiply  IV'  by  the  vtduo  of  Rayloiah  .sound  velocity  calculated  for  the  iivorii,!2,o  inititd  ptiramotins 


^..A  i.i;i  vtiUriffifil^ 


Part  ;6.  Test  cases,  resultis  auci  anscussioiis 


The -  t&sts.ofi  tile  ,the()re,tical  aioilei.  (levclopcd  winr  pciibniu'd  for  the*  rases  of  detouatiou 
propagatioii  in  iiydrogeii-oxygeii  luixtmi,*  in  a.  cylindrical  container  witli  tlnii  walls.  The 
par^neter.s-of  tlie  cpntahier  aud.tlu?  inixtun*  wcri,*.  <’ho.si*n  coire.spoiiding  to  the  experinients 
describeddn.  (4.1)..  Since  t lie  initi<d.,pre.sSuro. in  the?  tank  is  not  so  higii  it  wa.s  considered  that 
at  the  initiid  instant,  (f<  =  0  )  the  slKillJs  elastie*.  U.nder  fliis  as.smnption  the  iinitial  condilioiis 
(  t;  =„() )  for.  the  sliell  are  -the  followiii!* : 


0  =  0:  (7«.= 


Vojja 

'Ik  ' 


Sh  = 


=  O';  «i  =  r„£„\  r  =  (•(,  +  ((>; 


E  =  /, 


1  -r  exp(’2jfl)(2rfl  -  1)  //)i): 


D={)\  /)  = /)(d-2i«). 

wheie  t’o  -  initial  radius  <if  the  shell;  (c  -  displaccinent. 

Taking  into  account  tlu;  restrictions  on  the  minimal  fragments'  area,  oiitained  in  !j4.3  wi* 
a.ssume  that  for  the  sinalle.s't  fragments  their  area  is 

Totally  A'  ==  20  groups  of  fragments  wc'ie  loguarded  in  tin*  modt‘1.  The  areas  were  etiual 
to  S;  =  .Sj-i  •  j  =  2, 3,. ..20.  The  properties  of  the  materials  of  the  .shell  were  the 

following:  ft  -2700  kg/nv‘:  /t  =  27  GPa;  <>  =  924.3  .1/kg*  K:  a..  =  G.72  •  10-’*K-';  ./o  =  0.29 
GPa;  Jm«*  O.GS  GPa:  .i  =  125:  b  =  O.OCGGPa."':  \  =  G.2 •  IQ-'K-';  n  =  0.1. 

The  data  on  chemical  properties  of  comhustihle  mixture  inside  the  tank,  energy  release, 
combustion  rate,  activation  energies  and  DDT  proce.ss  is  shown  in  the  table  Fig.6.1.  The  data 
corresponds  to  the  hydrogen-oxygen  mixture  similar  to  tluit  used  in  tlu*  experiment  ESOC-2 

M;. 

Initial  conditions;  geometry  of  a  fuel  tank  and  physical  propert,ii*s  of  material  are  shown  in 
the  table  Fig.G.2.  Recalculated  propi'rties  of  mixture:  initial  and  final  eijuilihrium  states  ■  are 
shown  in  the  Fig.G.3.  The  key  parameters  of  detonation  wave  propagation  and  its  reflection 
from  the  walls  both  for  cylindrii'al  and  .si)lM*rical  tanks  for  the  giv<*n  initial  and  boundary 
conditions  are  shown  in  tlu!  Fig.G.4. 

The  pressure  -  time  diagram  of  wall  lotiding  is  shown  in  the  Fig.G.5.  It  should  be  mentioned 
separately  that  initial  time  (t  =  0)  on  the  diagram  for  the  case  of  {h*tonation  corresponds  to 
the  moment  of  the  onset  of  detonation  and  not  to  the  time  of  ignition  as  it  can  be  easily  seen 
from  the  table  (Fig.6.4). 

The  key  parameters  of  the  shell under  the  influenci*  of  the  loading  are  shown  in  the  table 
(Fig.6.6).  It  is  seen  that  both  for  ludindrical  and  spherical  shells  the  fragmentation  occurs 
under  such  condititnis. 

The  Fig.G.7  shows  the  diagram  of  tlu*  total  niuulu-r  of  fragments  of  difi’erent  mass  originating 
as  a  result  of  breakup  of  a  cylindrical  and  spherical  shell.  The  total  numlier  of  fragments  for 
cylindrical  case  is  larger  than  that  tor  a  spherical  one  due  to  the  jiarticular  lengthy^radius 
ratio  tor  the  cylinder  (see  Fig.G.2)  that  gives  a  largt*i  surface  in  comparison  with  a  sphere 
of  the  same  radius.  The  experinu*ntal  result.s  on  distrilmtions  of  fragments’  number  v(>rsus 
mass  are  shown  in  the  Fig.G.S  (exp(*rimental  data  [4.1]).  Comptirison  of  the  re.sults  shows 
that  for  the  part  of  spectrum  containing  large  fragnu*nts  coincidence  is  rather  fine.  A  slight 
difference  for  small  fragments  (tlu*ory  giv(*s  a  largi*r  number  than  the  experiment)  can  lu* 
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expitiiiiecl  1a'  tlu*  fact  that.iu  (-(iiirM*  «»f,  liu-t^xpcnimcut.N  not  all  thi?  iiagiacuts  were  collected. 
It  is  mehticnie*!  in  t4.l|  that  for  the  test  case  ES(’)C-’2  i.3‘X  , of  the  tottil  antss  was  lost  that 
equals  aljout  ,20  grams.  The  existing  (liffer«*nc<>  of  ma.ss  for  small  ficigiaeats  in  theoretical  and 
experimeatailfestdts  is  iimtdi  h\ss:  thiar  valtiei 

The  pl()ts  ()f  final  yelocityof  fragments  Versus  nnuss  for  diiiennit  angles  of  ofientatibh  oi 
fragments  are  givi-n  iii  the  Fig.G.9.  The  coincidence  of  fimil  vdocities  ;md;hreakup  velocities 
of  ffaghieiits  (Fig.G.G)  with  the  e.xpefihu*nEal  dtita  (4.1)  is  al.sts  fiii(‘. 

The  Fig.GTO  shows  the  total-monunittuu  dititrihiitioh  of  fnigments  versus  mass  for  fin- 
dric^  juicl.sphfefictil  ca.se.sTof  different  aiigles  oforieiitatioii. 

To  demonstrate  the  infiiveiice  of  the  type  of  tin*  process  of  energy  reletwe  iaside  the  tank 
the  next  .setof  figures  illustrtites  the  lin'tikup  of  the  same  fuel  . tank  in  cases  the  pure  cbmhus- 
tioh  ( deflagration j  t tikes  place  insitle  the  tank  {Figs.G.ll-G.lG)  and  deflagration  to  detonation 
trwisitibn  takes  place  near  the  wtills  of  the  ttink  (Figs.G;i7-G.22).  For  the  last  ca.se  aii  over¬ 
driven  (stfong).  detonation- wave  reflect.s,  from,  the  Willi  and  pressure  after  the  reflected  wave 
is  higher  (Figs.G.17.  GTS)  than  in  case  ofaiormal  detonation  (Figs.G.4i  G.5).  The  wall  loading 
for  i)DT  proce.ss  differs  for  cylindrical  iin< I  spherical  ciises  (Fig.G.lS).  .4s  a  re.sult  a  larger 
amount  of  elastic  energy  is  .stockpiled  liy  the  time  of  hrciikup  (Fig.G.lO)  iuid  the  number  oi 
fragments  increases  (see  Fig.Gi20).  \iitximum  number  of  fragments  is  moved  ill  the  direc¬ 
tion  of  smaller  masses  (Fig.G.20;  both  for  ludindricid  iind  spherictil  ctises.  The  final  velocity 
of  fragments  (Fig.G.21)  is  higher  th.in  in  ea.se  of  norniiil  detoniition  and  it  doesn't  depend 
practically  on  the  orientation  o  tt  fragment.  Total  monumtum  of  fragments  also  tunis  to  be 


liiglier  (Fig.G.22)i 

For  the  case  of  deflagration  wave  propagating  inside  the  tank  the  final  equilibrium  pressure 
in  the  tarik  is  the  same  as  for  detonatit)n  but  it  nei'ds  a  longin'  period  to  reach  tliis  pressure 
(Fig.6.12)i  The  re.sults  show  tlnit  by  the  time  of  breakup  elastic  energy  stockpiled  by  the 
wall  is  nearly  the  same  tis  that  for  the  .case,  of  normal  dett>nation  (Fig.G.13).  That  leads 
to  practically  the  same  fragmentation  pictun*  (Fig.C.14).  B\it  final  velocities  of  fragments 
(Fig.G.lS)  are  le.ss  than  in  ca.se  of  normtil  detonation.  Orientation  angle  of  fragment  in 
the  expanding  stream  of  reaction  products  plays  a  more  important  role  for  this  case.  Total 
momentum  of  fragntents  after  the  breakup  is  le.ss  than  in  case  of  normal  detonation  (Fig.C.lG). 
The  re.sults  show  that  the  angle  of  orientation  of  fragments  in  the  expanding  reaction  products 
plays  the  most  important  role  for  the  ca.sc  of  slow  loading  (di'fiagration)  (Fig.G.lG),  plays 
.smaller  role  in  case  of  normal  detonation  (Fig.G.lO)  and  practically  doesn't  influence  the 
results  for  the  ctcse  of  overtlriven  detonation  reflection  in  DDT  proce.ss  (p.G.21). 

The  results  of  numerical  modelling  show  thtit  the  Itrcidoii)  process,  the  number  and  mass 
distribution  of  fragments,  tmd  their  velocities  difler  greatly  depending  on  the  rate  of  energy 
relea.se  inside  the  tiuik  «md  this  dependence  is  not  monotonous.  For  the  lowest  (deflagration) 
and  highest  (detonation)  rtites  of  energy  releii.se  the  numl;er  of  fragments  is  less  than  for  the 
medium  rate  of  energy  release:  i.e.  deflagration  to  detonation  transition.  Maximal  fragments 
velocities  and  the  time  before  the  breiikup  for  those  three  c.i.ses  iue  illustriitetl  ])y  the  following 
table  for  cylindrictil  symmetry. 


Parameter 

Deflagration 

Tninsitjon 

proce.ss 

Detoniition 

Itreakup 
time  (ms) 

7017 

4272 

1029 

maximal 
velocity  (m/s) 

900 

2500 

1000 

Tliis  ()i  hfc'iikup  ('Itantctciisrirs  on  rlu'  nihf  of  tin*  (>ii(T"y  rcii'as(^  iuhidi'  thc-fucl 

tiuilC'  isstptv-lart^d  f<)  h(vu(>»lorttHl,  as  it.  was  usually  doiic  Ixd'oic. 

The  next  fit’iies  of  uiuiieni’al  c'xpoiimcuts  was dovorcd  ro  rcstiu!;;  the  so-ctdled  ">iinilarit.y 
pafmnetef>i"  for  Ijfd'aUup  (if'tikd  tank  (4.1).  The  .shiiiiarity  paran’ieter  was  stipposed  to  lie 
ihVo/M  T  4  =oousU  It  W!is  .supjkj.sc(l'  iu  (4.1)  rjiat  for  the  constant  walues  of  siniilarity 
paii’ametefs  the  iiniiiliH’  (jf  frasiuvnr's  phd  flioir  fiiiai  volocitic's  ichiainod  constant  and  masses 
(iflarge  ti%^'ment.s  ti;r('w  \ip  proportionally  to  tin*  »r(jwrh of  the  volmhe  of  a.slndl’s  matmiai. 

Figs.  G.23<  Gi24  present  tin*  resitl'tsof  calcuiatipns  for  the  cast*  similar  to  ESOC  2  1  ait  with 
pfoirpHiduully  iht'rkiusefl  radial  ditueiisious of  tin*  shell  and  thickness:  co  =.  130  sm:  h  =  2 
iniii.  (hi  ca.se  ES.QC  2  we  had-  /'»  =  30  .sni:  //  =  Om  mhn).  This  increase  of  dihiensionsof 
the  tank  iire.serves  the  .siinilhrity  parameter  .4  coiistant.  It  can  lie  seen-frt’hi  Fig.G.23  that 
the  iuimbeistif 'hYigments  ,c(>jitrary  to  predictions.  Tlie  incretiseof  hhal 

velocity  is  more  than  20%  (.see  Fig.6.24).  StiiaiTfraginents  have more  strong  flependehce  of 
velocity  on  the  ofieiitatidHv  angle. 

The  mass  of  large  fraginentS  diit*  to  the  i)h'dictions  was  tti  increa.st*  17  times,  hut  it  increased 
only  84=9  tihies  {conipiire  'Figs;G.7  atid  G.231; 

Thtis  the  resttlts  slitiw  that  there  cimitot  i'xi.st  siniilarity  paraineter.s  lor  .s\ich  si  complicated 
phenomena  as  hreiikuj) of  a  fuel  tank  aiid  .scjiling  of  model  ranks  is  a  very  ditticnlt  procedure. 
To  apply  the  n'Snlts  of  model  experinh'nts  to  riu*  real  fuel  taidc  it  is  iiecessary  to  make  use  of 
general  theoreticiil  model  of  the  type  w*  <l«*scril)i*d  iind  not  simple  sc.aling  parameters. 
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.-  A  'thepretidU  au<>44  t4V‘*  AuA -  tJUjlt  pt'  ;,;i)ii(’P  vciu<'i‘>  n^sult  of  (’luyU.i‘"'l 

expjpsiou.  is=  'vprkccL  put.  ..A  (\c)uipptpr  f<n^Inllnpll■i!.■<^l  .'jimulutiou of  tho  pro(!oss  is 

icpiistructecls  ^iie-4dsd(5asesiujit.iii-  tlie.exppniu<i 

The  results  of  nuinencal  iupflellih'^  slicny  tluit  r,li(.‘  hicakup  process.:  the  uvuuber.  iiiass 
clistriljutibasuaiu.iFdjIP'-'l'^id'’  (K.-poufliiis  on  xlie  race  pf  eiicrgy  release 

jusicle  theiaiik  for  cpU!jta,iit  ,t()tal  TXT  e<piiwal(,'uf  of  the  e.'cploshui.  The  lowest  rate  ot  energy 
.release  js  upriuahdefli^Sihtipu  (()r  coiuhustipu)..,the  iiu'diuiu  rate  -  ilpHagratibir  tp  fletouation 
tfaffsitioh,  theihighefiit  rate  -  <letonatioii.  Tt.is  prov.erl  theoretically  that  the  fiepeiuleiice  ot  the 
mailt  l)feakup,,paranieteis  on  tlai  rate  pleiiergy  relcase  is  nonlinear  ancl  it  is  not  nioiuitonous. 
■FIvus  it  is  necessary  to.  tajw  it  into  ac.couut  and  us<*  tliu  wirkecl  out  numerical  model  to 
calciilatedhei.irealdip  fpalm^^ 

it- is  proved  that  there cannot  ea.yist  .any  unitpu*  siniilarity  critericai  for  the  prblilem  ot 
breakup*<)f  n  fuel  tank  as  a  result  (nicheniicui  e.Kplpsuai.  The  breakup  .phenomenon  is  a 
combination  of  complicated  prp<*essesipf  .(lifteienr.  pliysicad  natvue.  That  is  why  there  exist  a 
lotof.diniensionless  .siniilarity  parame,ier.s  audit  is  impaissilile  to  satisfy  till  of  thein  without 
c6mihg.to.il  c.ontradictiiiii.  Thus  .sniling of  mixlel  <'xperimeuts  is  impossible  withoiit  makinii 
iise  of  a  profound  theoretietil  iind  numerical  mod«i  of  the  type  descrilietl  in  this  report. 

Some  of  the  important  stiiges  of  the  process  cati  be  the  topics  of  further  investigations. 
More  complicated  geometries  of  the  shell  will  leiul  to  new  problems.  The  presence  of  weak 
zones  in  the  shells  of  variable  thicknc.ss  can  lead  to  a  eomphiely  different  seenafio  of  breakup 
in  ca.ses  of  slow  and  intense  loiulings.  The  lireiikup  seemirio  for  different  combinations  of  the 
variations  of  shell’s  thickiie.ss  (shell  uommiformities)  and  rates  of  energy  release  is  also  to  be 
investigated; 
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Editing  reaction  s  data,  <ESC>  to  confirn 

- Input  data  on  the  chenical  reaction  stoichionetry - 

Brutto  reaction  elenents:  XI  2  ♦  XI  0  ♦ 

XI  -  fuel,  X2  -  oxidizer,  XZ  1  ♦  X2  0  ♦ 

X3  -  reaction  produt,ts,  X3  0  »  ==>  X3  2  ♦ 

Xt  -  neutral  conponent.  X3  0  X3  0 

- Input  data  on  the  reagent's  properties - 

tiolar  nass  (g/nol): 

XI  2.00  XZ  32.00  X3  18.00  X4  28.00 

Initial  uolunctric  concentrations: 

XI  .6667  X2  .3333  X3  .0000  X4  .0000 

Initial  adiabatic  ratio: 

XI  1.40  X2  1,40  X3  1.33  X4  1.40 

Final  adiabatic  ratio: 

XI  1.33  X2  1.30  ■  X3  1.20  X4  1.32 

- Input  data  on  the  reaction  energetics  and  ignitiuity - 

Reaction  heat  rate,  kJ/nol  Octluation  energy,  kJ/nol 

118.000  70.000 

Predetonational  length,  cn  Deflagration  speed  rate,  cn/s 

12.000  1.000 

Unsteady  detonatlonal  length,  cn  Detonation's  fornatlon  length,  cn 


XI  1.33 


X4  1.40 


X4  1.32 
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Editing  conditions  and  geonetry  data,  <ESC>  to  confirn 

- input  data  on  initial  internal  conditions - 

Initial  pressure  (bar)  4.180  tenperature  (K)  Z30.0 

- Input  data  on  geonetry - 

Fuel  tank  radius  (cn)  Z5.e  length  (cn)  1Z5.0 

- Input  data  on  the  nenbrane  uall  geonetry  &  thernodynanics - 

Thickness  Density  Heat  capacity  Uolune  extendibility 

«n:  .58  g/cn*3:  2.700  J/(kg*»K):  924.3  (1/HK)  67.20 

- Input  data  on  the  stresses  nodules  d  uiscosity - 

Shifting  nodule  Elasticity  Unit  tiax.  elasticity  Unit 

6Pa:  27:O00  GPa:  .298  GPa:  .680 

Nenbrane  dynanical  uiscosity  (kN»s/n)  108.0 

- Input  data  on  Steinberg  -  Gu inane  nodel  constants - 

Beta  125.0  b  (l/6Pa)  .065  h_(l/kK)  .620  n  .100 

- Input  data  on  destruction  criteriun - 

Haxinal  dissipation  (kJ/kg)  30.808 

Energy  per  breach  square  unit  (kJ/n^Z)  108.000 

Fragnentation:  charact. square  (cn^2)  4.000  expon.paraneter  .50 
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Internal-  process  results 
■Concentrations  of  reagents 


XI 

X2 

X3 

X4 

Uolunetric  initial  concentration 

.'66670 

.33330 

.00000 

.00000 

ilass  initial  concentrations 

.11112 

.88887 

.00000 

.00000 

Uolunetric  final  concentrations 

.00014 

.00000 

.99985 

.00000 

Hass  final  concentrations 

.0000,1 

.00000 

.99998 

.00000 

- Holar  nasscs  of  reagent  nixtures,  g/noi - 

Initial:  11.9  Final:  17.9 

- ^Specific  heats  of  nixture  at  constant  uolunci  J/(kg*K) - 

Initial:  1731.3  _  Final:  2308.5 

- Adiabatic  ratli  of  nixture - 

Initial:  1.400  Final:  1.200 

- Reaction’s  heat  release - 

Release  per  nol  (kJ/nol):  117.988  Release  per  nass  (kJ/kg):  9833.169 

- Equilibrlun  final  conditions - 

Pressure  (bar):  52.673  Tenperature  (K):  4432.0 

- Type  of  the  process  inside  the  fuel  tank - 

Nornal  detonation 
Press  any  key  to  continue 
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Internal  process  results  (continued) 
Tine-pressure  characteristics  of  the  process 


Pressure  (bar): 

Cylinder 

Sphere 

Initial 

1.100 

1.100 

Before  shock 

1.100 

1.100 

After  shock 

191.385 

191.385 

After  reaction  zone 

93.061 

93.061 

After  shock’s  reflection 

1376.750 

1376.750 

After  reaction  zone's  reflection 

2Z9.173 

229.173 

Final 

52.673 

52.673 

Detonation's  uauc  velocity  (kn/s)^ 

3.015 

3.015 

Tine  intervals  (ns); 

Ignition  —  detonation  origin 

1028.829 

669.182 

Detonation  origin  —  reflection 

.013 

.013 

Reaction  zone  near  the  uall  duration 

.009 

.009 

Relaxation  tine 

.082 

.082 

Active  pressure  Increase  duration  .09Z 

Press  any  key  to  continue 

.092 

Tank  uall  breakup  by  internal  pressure  increase  results 
- Loading  on  the  uall - 


Inpulse  per  square  (bar»s) 

Cylindrical 

.011871 

Spherical 

.011871 

Dissipation  per  nass  (kJ/kg) 

36.5492 

36.5492 

Elastic  energy  per  nass  (kJ/kg) 

78.7556 

78.7556 

Final  density  (g/cn''3) 

1.9131 

1.9131 

Final  uelocity  (n/s) 

973.4301 

973.4301 

Loading  tine  (ns): 

.2001 

.2001 

Uall  destruction  tine  (ns): 

.1232 

.1232 

Fragnentation  occurs: 

Ves 

Tes 

Press  any  key  to  continue 
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Internal  process  results  (continued) 
■Time-pressure  characteristics  of  the  process 


Pressure  (bar); 

Initial 

Before  shock 

After  shock 

After  reaction  zone 

After  shack’s  reflection 

After  reaction  zone's  reflection 

Final 

Detonation’s  uaue  uelocity  (kn/s) 

Tine  interuals  (ns): 

Deflagration  process  duration 


Cylinder 

Sphere 

4.100 

4.100 

4.100 

4.  ICO 

4.100 

4.100 

4.100 

4.100 

4.100 

4.100 

4.100 

4.100 

52.673 

52.673 

.000 

.000 

7017.515  5296.636 


Press  any  key  to  continue 


Tank  wall  breakup  by  internal  pressure  increase  results 
- - Loading  on  the  uall - 


Inpulse  per  square  (bar«s) 

Cylindrical 

Spherical 

.010069 

.010069 

Dissipation  per  nass  (kj/kg) 

37.4761 

37.4761 

Elastic  energy  per  nass  (kJ/kg) 

78.5099 

78.5099 

Final  density  (g/c«''3) 

1.7791 

1.7791 

Final  uelocity  (n/s) 

877.4780 

877.4780 

Loading  tine  (ns); 

.2871 

.2871 

Uall  destruction  tine  (ns): 

.1911 

.1911 

Fragnentatlon  occurs: 

Ves 

Ves 

Prc.ss  any  key  to  continue 


Fragnents  final  ueiocity 

y  (kn/s)<  cylindrical  U  (kn/s)i  spherical 


Orientation 


Fragnents  final  total  impulse 

M"y  (kg*<fi/s)i  cylindrical  (kg*ti''s).  spherical 


i.n 


60 


«,V' 


Tank  uall  breakup  by  internal  pressure  increase  results 


- Load  mg 

on  the  uall - 

Cylindrical 

Spherical 

Inpulse  per  square  (bar«s) 

.031508 

.026749 

Dissipation  per  nass  (kJ/kg) 

33.1112 

33.5926 

Elastic  energy  per  nass  (kJ/kg) 

94.7972 

89.5447 

Final  density  (g/cn*3) 

2.1292 

2.1175 

Final  uelocity  (n/s) 

2580.6010 

2172.9766 

Loading  tine  (ns): 

.0590 

.0666 

yall  destruction  tine  (ns): 

.0337 

.0366 

Fragnentatlon  occurs:. 

Tes 

Ves 

Press  any  key  to  continue 
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Nunbcr  of  particles  diagramnes 
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Fragments  final  uelocity 

U  (km/s).  spherical 


U  (kn/s)i  cylindrical 


Fragff 


M"V  (kg**n/'s)i  cylindrical 
0-) - , - , - , - 


;  final  total  inpulse 

H“U  (kg*»ti/s),  spherical 
1000  J - ^ - - 


Orientation 


Appendix.  L'ser’s  instruction. 

The  program  calcukucM  picfliction  of  fla“uu•lu^  ntuiii)i  i.  !ua>.s-imiiibei  (iistnbuiiou.  hiial 
velocities  and  moineiituiu  (lepeu<iiu!*  on  i()mi)nstioii  pioccss  lu  the  closed  vessel  of  sim¬ 
ple  shape  -  eylindiic.il  with  hemispluTical  ed“es.  The  theoreticjil  pait  of  calculaiion  is 
described  i  i  the  repoit. 

The  program  retpiires  IBM  .A.T  or  eompafii)le  3SG/S7  or  4SG  witii  \’GA/SVGA  monitor. 
MSDOS  ver.3.30  or  higher  tnid  <)40K  memory. 

Minimai  .set  of  files: 

FTB.EXE  -  program  itself 

14F0NT  -  font  for  “laphie  applictitioiis.  In  ctise  ot  alssceiiee  of  this  tile  the  program  in 
graphic  modes  will  use  deftuilt  text  chartieter  font  Init  with  less  ticcnracy. 

Extended  set  of  files: 

*.INI  -  user-named  files  eontaining  iiiititil  data  infornuitioii. 

*.RES  -  user-named  tiles  eontaining  re.sults  of  ealcuhitions. 


Proee.ss  of  ealeulations  deseiiptioii. 

0.  PRELIMINARY  NOTES. 

The  program  is  controlled  only  by  keylxiard  with  no  use  of  a  mouse  To  print  ."esults  some 
extei'nal  resident  program  is  iietxled  (stu'h  ;i  piogitun  mostly  uses  <  Print  .sereen>  key  to 
be  activated).  If  you  want  to  ttdte  some  piintetl  copies  from  the  plots  yoti  must  run  such 
a  program  before. 

1.  START. 

To  start  the  progrtun.  one  should  run  FTB.EXE.  .\ft('r  tht'  first  screen  appears  one  ctui 
press  <ESC>  to  exit  immediately  or  any  other  key  to  enter  tlie  INTTI.4L  .MENU.  If  the  file 
14F0NT  doesn't  exist  or  tin*  system  ciinnot  access  it  or  it  theiv  is  l(>ss  available  memory 
than  needed,  then  tlie  waniing  inessatit'  appetirs  on  scieeii: 

-  ’’Cannot  acce.ss  file  or  allocate  memory  for  graphics  tout  diitiT 
Then  the  system  anyhow  changes  to  INTTI.\L  MENU. 

2.  INITIAL  .MENU. 

This  menu  contains  3  items: 

DATA  FROM  FILE.  INPUT  DEFAULT  D.4TA  ;uid  EXIT  TO  DOS.  Ntivigation  Ix'tween 
items  is  performed  with  <Right>  and  <Left>  keys.  To  choost*  an  item  one  should  press 
<Enter>  key.  The  <Esc>  key  heie  is  in:ivaila!)le.  The  piimarilv  active  item  alwavs  is 
INPUT  DEFAULT  D.\TA. 

2.1.  DATA  FRO.M  FILE. 

The  system  asks  for  the  fih*-  name  ( withoiu  extension )  wiu'ie  to  use  initial  data  from.  Then 
it  tries  to  open  this  file  with  .INI  extensioirtmd  to  trad  initial  (kita  from  it.  If  this  file 
doe.su't  exist  or  the  access  is  denied  then  the  system  says  "Ciitinot  acce.ss  hie”  and  remains 
in  the  INTTLA.L  MENU.  Else*  it  reads  the  initial  data  and  changes  to  INITL4L  DAT.A. 


EDIT.  Note  :hut  the  tilciiamc' then  tluu^l  to  ])<■  acrivc 

2.2.  INPUT  DEFAULT  DATA?  '* 

The  .system  just  simply  ('h:iui;,es  to  INITIAL  DATA  EDIT.  Tlu*'  active  initial  data  idefatilt 
or  previously  read  from  a  file)  i^ets  chtmee  to  he  edited  for  new  ctdcuhition. 

2.3.  EXIT  TO  DOS.  (See  FINISH). 

3.  INITIAL  D.\TA  EDIT. 

.A.t  first  the  system  opens  :i  screen  to  edit  the' followinu,  list  of  dtita  (first  input  screen): 

-  data  on  the  chemietd  reaction  f.toichiometry:  stoicliioiuetrir  coefficients  of  the  chemical 
reaction  in  the  vessel,  initial  tmd  final,  for  the  set  of  u,enertdized  components  -  ftiel.  oxi¬ 
dizer,  reaction  product  and  netitrah 

-  data  on  the  reagent's  properties:  moltu  imisses.  mititd  volumetric  concentrations,  initial 
adiabatic  ratio  and  final  adiabatic  ratio: 

-  data  on  the  reaction  energetics  tmd  ignitivity:  leaction  heat  release,  tictivation  energy, 
predetonational  length,  defiagration  spee<l.  unsteady  detonation  length,  detonation  forma¬ 
tion  length. 

Purposes  for  this  data  are  described  in  the  leport.  To  choose  a  field  to  edit  one  can  use 
<Left>  .  <Right>.  '-.Up>.  <Down>,  •;,Home>  and  <Encl>  keys:  note  that  <Left>  and 
<Up>  do  the  stune  its  like  as  <Riglit>  and  <Down>.  To  Ijigin  editing  a  field  press 
<Enter>  or  any  other  (no  mivigation's)  law.  To  confirm  (sliting  and  exit  press  <Esc> 
Then  if  some  limitations  tire  violated.  :i  w.uning  message  tippears  on  the  screen,  ami 
editing  continues.  When  ti  Held  is  in  the  process  of  editing.  <Left>  tmd  <Right>  moves 
the  cursor.  <Del>  and  <D!ickspace>  performes  symbols  deletion.  <Enter>  confirms  and 
<Esc>  lueaks  editing. 

Limitations  for  the  vidues  edited  on  the  first  input  screen  produce  the  following  warning 
messages  in  case  of  their  violation: 

-  ’’Invalid  data  on  reaction  coefficients" 

(stoichiometric  coefficients  must  be  non-negtitivei: 

-  "No  fuel  on  reaction’s  inimt" 

(fviel  stoichiometric  coefficient  before  the  reai-tion  mu.'-t  be  positive i: 

-  "No  inoducts  on  reaction's  output  ’ 

(products  stoichiometric  coefficient  aftei  the  letiction  must  Ite  positive!: 

-  "Netral  component  can't  change  in  reliction " 

(neutral  component  stoichiometric  coefficients  befoie  and  after  the  reaction  must  be  the 
same); 

-  "Invalid  molar  mass  of  a  component  particijiating  the  retiction" 

(molar  ma.ss  of  each  coniiiouent  with  positive  stoichiometric  coefficient  must  lie  po.sitive): 

-  "Negative  or  zero  initial  concentrtition  of  an  inimt  component 

(initial  concentration  of  each  coiniionent  with  positive  stoichiometric  coefficient  before  the 
reaction  must  be  positive!: 

-  ".4diabatic  ratio  of  participating  comitonenfiiiusi  l»e  from  1  to  1.G7". 

-  ’’Molar  mass  of  components  must  lie  the  same  both  on  input  tmd  outimt' 

(total  molar  ma.ss  of  components  must  be  the  same  both  before  tmd  after  the  reaction!: 
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-  "Sum  of  voi  .uietiic  coiK'omiatious  ou  input  uium  in-  I  . 

(reaction's  heat  must  lie  positive): 

-  "Negative  or  zero  activation  energy  " 

(activation  energy  must  he  positive): 

-  "Negative  or  zero  predetonatiomil  length  " 

(predetonational  length  must  he  positive); 

-  "Negative  or  zero  deflagration  speed  rtite  " 

( deflagrational  speed  must  lie  positive): 

-  "Invalid  detonation's  formation  length" 

(detonation's  formation  length  must  he  positive  lint  not  exceeding  piedetonational  length); 

-  "Invzdid  overdriven  detonation's  fonmition  length " 

(overdriven  iletonatioii’s  formation  length  must  lie  positive  hut  not  exceeding  detonation  s 
formation  length). 

Then  the  system  opens  a  screen  to  edit  the  following  list  of  dtita  (sc'cond  input  screen): 

-  data  on  initial  internal  conditions:  [iiessure  and  temjierature: 

-  data  on  geometry:  fuel  t;mk  radius  and  length.  Note  tlnit  "length"  here  means  ttital 
length  of  the  vessel  including  hemispherical  eclges: 

-  data  on  tla*  memhriuie  wtill  geonietiy  ami  thermodynamics;  thickness,  density,  heat  ca¬ 
pacity  and  volume  extendiliility  (with  heating); 

-  data  on  the  stresses  modules  and  viscosity:  shifting  module,  elasticity  limit  (initial), 
elasticity  limit  (maximal),  dynamical  viscosity  of  the  wall  nuitmial: 

-  data  on  Steinherg  -  Guiiiane  model  cousttuits:  -I  values  deiiending  on  the  wall  imuerinl: 

-  data  on  destruction  criterium:  maximal  dissipation,  energy  per  lireiich  stpuire  unit,  char- 
acteristictd  stpuire  of  fiagments  tnid  exiionential  iiarameter  in  the'  distrihution  law. 

All  rules  of  eiliting  are  the  same  as  for  the  first  input  screen. 

Limitations  for  the  second  input  scteeii  pioduce  the  following  warning  messiiges  in  case  of 
their  violation; 

-  "Invalid  initial  jnessme  value" 

(initial  pressure  must  lie  positive): 

-  "Invalid  initial  temperatuie  vtilue" 

(initial  temperatuie  must  l>e  not  helow  2()K  and  aliove  .j(JOl)K); 

-  "’Fuel  tank  length  must  lie  etpial  oi  exceed  two  radii" 

(to  perform  the  vessel  with  :t  middle  cylindrical  part  and  two  hennspheiicjil  edges); 

-  "Tank  wall  thickness  less  oi  eipuil  to  zero" 

(thickne.ss  of  the  tank  wail  must  lie  iiositive); 

-  "Wall  density  must  oe  iiositive": 

-  "Wall  heat  cainicity  must  lie  positivi*"; 

-  "Wall  thermal  volume' extendiljility  must  lie  positive". 

-  "Wall  shifting  module  must  lie  positive": 

-  "Wall  elasticity  limit  must  lie  iiositive’"-; 

-  "Maximal  wall  elasticity  limit  must  exceed  the  iniiitd  one"; 

•  "’Wall  dynamical  viscosity  must  lx-  positive"; 

-  "Each  of  the  Steiifneig  Guinane  model  constants  must  lie  oositive  : 
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-  "Maximal  di-ssipatiou  nnihr  he*  poMtivi'": 

—  ■ . .  . . 

-  '’Charactcristical  scjuaie  of  fraiimoutation  umst  he  posinvo": 

-  '’Negative  exponential  parameter  m  fitigmentation  distrilaition " 

(it  must  be  (lositive  or  zero). 

After  confirming  editing  of  the  second  input  scieen  the'  system  clumges  to  the  MAIN 

MENU. 

4.  MAIN  MENU. 

The  MAIN  MENU  of  the  prognun  consi.sts  of  9  items.,  change  between  them  is  performed 
with  <Left>.  <Right>.  <Up>.  ■:,Down>.  <H()me>  and  <End>  keys.  Item  is  chosen 
with  <Enter>  key.  <Esc>  key  is  iiuivailable.  When  one  enters  the  MAIN  MENU  from 
the  INITIAL  DAT.\  EDIT,  the  inititi!  itc'in  is  C.\LCULATE.  Items  of  the  main  menu  are: 
CALCULATE,  NEW  CALCULATION.  EXIT  TO  DOS. 

SHOW  RESULTS.  S.WE  RESULTS.  READ  RESULTS. 

SHOW  INITIAL  D.\TA.  S.WE  INITIAL  D.\TA.  SERX’ICE. 

4.1.  CALCULATE. 

The  system  calcuhites  results  if  the'  Hag  "c.dculated"  is  off.  This  Hag  is  off  on  enter  from 
INITIAL  D.4T.\  EDIT  <iud  turns  on  td'tei  calculation  oi  RE.4D  RESULTS  (in  ca.se  of 
succesful  reading).  When  this  Hag  is  on.  choosing  of  this  itiun  letids  to  warning  message: 
"Results  for  the  active  initial  data  luive  lieen  ahetidy  olitained".  ;ind  the  choice  pseudo¬ 
cursor  in  M.\IN  MENU  turns  to  SHOW  RESULTS.  Cidculation  consists  of  three  ste  ,“s: 
internal  proce.ss  (deflagrtition.  detonation,  ddt-process  or  no  process  at  all],  wal)  '  '  ..ing 
and  fragmentation.  Hnal  velocities  calculation.  If  the'  H.ig  "output  during  calculations"  is 
on  (it  can  be  switched  on/off  in  SERVICE),  then  after  each  stage  temporary  results  will 
be  displiiyed  (all  these  n'stilts  can  id.so  be  displayed  via  SHOW  RESULTS).  By  default  the' 
Hag  "output  (hiring  ctdculations"  is  oH'. 

4.2.  NEW  CALCUL.\TION. 

The  system  turns  tigtiin  to  LNTTL-\L  MENU.  Note  thtit  the  Hag  "calculati'd”  turns  off.  but 
the  active  inicitd  ihita  is  conscnvisl  as  well  as  tin*  .active  nleniime. 

4.3.  EXIT  TO  DOS.  (See  FINISH). 

4.4.  SHOW  RESULTS. 

If  the  H.ag  "internal  details  output"  (--ee  SER^’1CE)  i^  on.  then  the  systiuu  shows  the 
following. 

First  screen  of  internal  lesults. 

-  Concentrations  of  retigents; 

-  Volumetric  initial  concentiaiions: 

-  Mass  initial  concentr;itions: 

-  Volumetric  Hnal  concentrations: 

-  Mass  Hnal  concentiations: 

-  Molar  mas.ses  of  leagent  nnxtuies  (initial  and  Hiudi, 


-  Reaction's  heat  release  (per  mol  and  pei  mass  uiiir  i; 

-  Equilibrium  feud  conditions  (pressure  and  temp<-iatiue ): 

-  Type  of  the  process  inside  the  fuel  tank  (no  piocess.  deHaiiiarion.  normal  detonation, 
ddt  in  pro,e,ress  or  ddt  oriainatin** ). 

Second  screen  of  internal  residts. 

For  spherical  tuul  cylindiic;'.!  i)arts  of  the'  v<'sstd  the  tollowma  is  shown; 

-  Pressure  (initial.  Indore  shock  wave,  altt'r  shock  'wave,  after  iea<'tion  zone,  after  shock  s 
reflection,  after  reaction  zone's  (‘flection  and  final i.  In  case  of  "no  piocess"  or  ''deflae;ra- 
tion”  only  initird  and  final  pre.ssure' values  have' resouaide  meanina. 

-  Detonation's  wave  velocity  (for  detonation  oi  dilr  proce.ss).  -  Time  intervals  (ij;,nition 

-  detonation  orioin.  reaction  zone  near  the  w;di  duration,  relaxation  time,  active  pres¬ 
sure  increase  duration).  Some  inii  ivals  have  ri'sonalile' meaiiina,  only  in  detonation  or  ddt 
processes. 

Third  st'reen  of  internal  process  shows  the  chart  of  piessure-time  evolution  near  the  wall 
(in  graphics  mode).  Data  for  cylindricaLcase  is  shown  m  white,  foi  spherical  -  in  retl.  If 
the  curves  (are  the  same  (i.e.  when  tioimal  detonation  oci  uisi  the  led  cnive  is  drawn  ovi'r 
the  white  one. 

If  the  flag  "meiiilnane  Inirst  details  output "  is  on  (see  .SER\’ICE)  then  the  system  shows 
the  following. 

First  screen  of  meiiiljraiie  burst  and  loading  outjuit. 

If  no  proce.ss  in  the  vessel,  then  the' systtnu  simply  ty]>es  “No  cheniiccd  retiction  in  the  ves¬ 
sel”,  else  for  spherical  and  cylindrical  pin  ts  of  the  vessel  the'  following  results  are  displaj'ed. 

-  Impulse  per  sipuue: 

-  Di.s.sipation  per  iiiiiss  unit: 

-  Elastic  ener,g\'  per  nutss  unit; 

-  Final  density: 

-  Final  velocity; 

-  Loading  time  (time  of  piessure  iqipliciition ): 

-  Wall  destruction  time  (in  ciise  of  bieakuii  it  is  usiudlv  less  then  the  Loading  time',  othei- 
wise  it  is  zero.  i.e'.  it  has  no  reasonable  meaning i. 

-  If  fragmentation  occurs  (yes  or  no). 

Note  that  the  first  four  values  iue  cidcuhited  until  the  bieaku])  occiiis  or  the  conditions 
show  that  it  will  not  occur:  the  5th  (final  velocity;  vtdue  in  I'tise  ot  Ineakup  lelates  <  the 
moment  when  the  huge  increase  of  piessuie  uetu  the  wall  will  lelax  and  the  fragments  then 
will  accelerate  mostly  liy  the  driig  forces. 

Second  screen  of  memlntiiie  burst  and  lotiding  output. 

Shows  ditigi'iimmes  of  mas.s-numliei  distrilmtion  of  fragments  foi  cylindiicid  case  (in  gray) 
and  spherical  ca'-e  (in  dark  led).  If  the ’diiigniniines  tor  both  ctises  titinsect  then  one  will 
see  them  m  pink.  Diagrannnes  are' shown  m  giajihii's  mode  in  logtirithmic  axes. 


If  the  fiau,  'final  vrlocirics  output "  is  on  tscc  SERMC'Ei  then  the  system  sliows  the  fol- 
.  . . . . . . .  . . 


First  screen  -  final  velocity  dependiin!,  on  iia>:,nient  s  mass  plots  foi  cyiindiictd  and  spherictd 
cases  are  shown.  Different  emves  lelate'  to  different  taii'iication  ;inu,ies  of  frai>,ments:  they 
are  shown  in  different  colors.  ;ind  the  corre-spondinu  .unties  tire  shown  in  the  ri2,ht  ptirt 
of  the  screen  eacii  in  its  own  color.  If  no  iiamnentarion  oecuis  for  some  ctise  (cyiindiictd. 
spherical  or  liotli  of  them)  then  "No  fra»mentation "  is  typed  instead  of  a  plot. 

Second  screen  -  final  momentum  of  frattments  dependiim,  on  mass  plots  for  cyiindiictd 
.and  .spherical  ca.ses  tue  shown.  Different  curves  lehue  lo  diffeieut  oiientation  ambles  of 
fragments  like  on  the  pievious  screen. 

Note  that  change  hetwt'en  screens  is  intide  liy  piessimi,  ti  key. 

4.5.  SAVE  RESTLTS. 

This  item  of  the  MAIN  MENU  works  when  "calculated  “  Hag  is  on  and  the  active  filename 
is  determined.  Then  the  lesults  are  saved  into  Hie  ntmietl  <Hlenanie>.i'es.  Otherwise  the 
following  nu'ssages  can  tippetir  on  screen: 

-  "'Results  on  the  tictive  initial  data  litiven't  heen  olitiuned " 

(the  flag  "calcuhited"  is  oH  k 

-  "Cannot  access  Hie.  write  initial  d.ita  at  Hist" 

(filename  is  undefined  oi  cannot  open  Hie  with  deHiied  name). 

4.6.  READ  RESULTS. 

This  item  of  tlie  M.4IN  MENU  works  when  the  active  Hlemune  is  determined.  Then  the 
results  .lie  retid  fiom  Hie  named  <Hlen<'inn‘>.ies.  Otheiwise  the  following  message  can 
appear  on  scrmi; 

-  "Cannot  acct*ss  Hie” 

(Hlename  is  undefined  or  eaniioi  open  Hie  with  defined  name). 

Note  that  succesfull  ptiss  tluough  this  item  of  the  .\I.\iN  MENU  sets  '"calculated”  flag  on. 

4.7.  SHOW  INITIAL  D.\TA. 

The  active  initial  data  is  shown  in  the  oidei  descuhed  ahove  m  INITI.\L  DATA  EDIT, 
but  editing  is  now  iiiavailalile. 

4.S.  SAVE  INITIAL  D.-\TA. 

The  system  asks  for  the  file  name  (without  e.xti'nsiou ).  Tlie  active  iilenanK*  is  default  if 
determined.  Then  it  tries  to  open  a  file  nameil  < Hlename mi  and  to  save  initial  data  into 
it.  If  opening  file  is  unsuceesfull  then  a  waiimig  messaiic  appears: 

-  "Cannot  access  Hie". 

Note  that  in  case  of  succesfull  writ»‘  to  Hie  th<-  active  Hlename  is  i  ie)deti'rtnined. 

4.9.  SERVICE. 

This  item  of  the  M.-VIN  ?\IENU  leads  to  sul  )-menu  I’oiisistinn  ol  H\’e  items; 

-  "Proceed  (Bypass)  interiiai  details  ouijutt '. 


-  "Procc(-“(l  iByp.iss)  iiu'iui)i;iiic  i)msr  (iriniK  onruiu'  . 

‘-~*T’TbV'w(TtByiTn.'<sTtiYr'fhrctoiwiT*v(‘wfT7nr'f'^--  •  , 

-  "Pi’ocwd  (Bypass)  ()\itpm  dmiii!*  call  ularions  . 

-  "Exit  to  lujiiii  menu". 

Each  of  the  first  four  items  si'rves  to  switch  ou/otf  the  l•olrespolldinu,  contro!  iiaii.  the 
.status  of  each  fltie  is  shown  in  the'  piompt  (i.e.  if  "Pioceed"  is  in  the  prompt  then  the 
flag  is  on).  The  purposes  of  the.se'  Htigs  can  lie  seen  tiixwe  in  CALCULATION  and  SHOW 
RESULTS.  One  can  exit  to  imiin  menu  using  the  nth  item  or  just  pressing  the  <,ESC> 
key.  Note  that  if  the  tluee  first  flags  aie  OFF  then  .SHOW  RESULTS  will  not  show  any 
results  fit  all. 

5.  FINISH. 

The  first  scieen  fippears  :ig:iin.  Lut  with  "P:es>  ;inv  key  to  exit"  on  it.  Then  aftei  <i  key  is 
pressed  the  program  exits  to  DOS. 
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